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The doublet character of the infra-red absorption line at 2.674 of solid potassium hydrogen 
fluoride is due to the “tunnel effect.” The proton has two potential minima between the 
fluorine nuclei in the FHF~ ion. The distance of a minimum to the center of the ion is 0.26A. 
The height of the potential hill is 33,400 calories. 


HE infra-red absorption line of solid potas- 

sium hydrogen fluoride at 2.674 has been 
found to have doublet character by Buswell, 
Maycock, and Rodebush.! The two lines are 
located at 2.6494 (3775.0 and 2.729u 
(3664.4 cm). Buswell et al. suggested that this 
doublet nature might be due to the ‘tunnel 
effect” on the supposition that the proton 
may oscillate between two potential minima 
located between the two fluorine nuclei and 
equidistant from a plane at right angles to the 
figure axis of the FHF- ion at its midpoint. 

A similar case in the ammonia molecule has 
been studied by several authors?~‘ with a variety 
of potential energy expressions. The values of the 
total energy given by Wall and Glockler have 
been applied to the present case, after having 
solved them for all possible values of the param- 
eter g. This quantity is given by 


where uw is (1/38+1/1.008) or 0.982, vo is the 
average fundamental frequency (3719.7 cm=), / 
is the distance from the center of the ion to the 


1A. M. Buswell, R. L. Maycock, and W. H. Rodebush, 
J. Phys. Chem. 8, 362 (1940). 

*F. Hund, Zeits. f. Physik 43, 805 (1927). 

*D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 

*See F. T. Wall and G. Glockler, J. Chem. Phys. 5, 314 
(1937) for further literature references. 


potential minimum, and h is Planck’s constant. 
For the purpose of calculation the following 
abbreviations have been used: 


exp 


B=4/(m)'g-exp 

C=exp (—@’); 

D=(1—29) exp (—’). 
Wos/hvo=3[1+(2¢?—A —B)/(1+C)]; (2) 
Woa/hvo= (2¢?—A)/(1—C)]; (3) 
Wrs/hvo=3(1+(2g?—A —Bq’)/(3—3D)]; (4) 
—2B+Bq’)/ 

(3+3D)]. (5) 

The calculated values are shown in Fig. 1. 
The doublet separation is given by 

Ahvo= Wra— Wist+ Woa—Wos (6) 

and is plotted as a function of g in Fig. 2. 

From this graph a value of g can be found from 

the experimental doublet separation Ahvo/hvo 


= (3775.0— 3664.4) /3719.7 =0.03. In the present 
case g is 2.7. From Eq. (1), which can be written 


1=5.805q/(uv0)?, (1a) 


the distance (1) from the center to the potential 
minimum is found to be 0.26 angstrom. 
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Fic. 1. Energy states (W) for double minimum problem. 
S=symmetric and A=antisymmetric levels; 
mental frequency ; g=coordinate. (Applicable for values of 
the parameter g greater than about 2.) 


In order for the whole consideration to be 
valid this distance must check reasonably well 
with the dimensions of the FHF- ion as obtained 
from other sources. The following information 
has been found on the size of this ion. Bozorth® 
studied the diffraction of x-rays from potassium 
hydrogen fluoride and found the F—F distance 
to be 2.25 angstroms. He could not obtain the 
H—F distance since protons usually do not 
scatter x-rays sufficiently to be observed. The 
distance from the center of the ion to a fluorine 
nucleus is then 1.13A. Another value may be 
obtained from the application of the theory of 
small vibrations to the triatomic FHF- ion. 
The line at 2.67 is considered to be the asym- 
metric frequency : 


v3=[(k.p)/m ]*=3719.7 


where k=force constant, m=mass of fluorine 
atom, M=massof hydrogen,and p= (2m+M)/m. 
The simple valency force field has been used.® 
The force constant calculated is 0.4 megadyne 
per cm. An application of Badger’s rule’ to the 
hydrides of the first row of the periodic table 
shows that the rule holds very well and permits 
the calculation of the H—F distance in FHF- 
from the force constant. The value*® obtained is 


5 R. M. Bozorth, J. Am. Chem. Soc. 45, 2128 (1923). 

®K. W. F. Kohlrausch, Der Smekal-Raman Effekt 
iy (Julius Springer, Berlin, 1938), p. 64. 
(1935) M. Badger, J. Chem. Phys. 2, 128 (1934); 3, 710 
8 Unpublished results. 
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Fic. 2. Doublet separation in tunnel effect. »»o=funda- 
mental frequency ; g=coordinate. (Applicable for values of 
the parameter g greater than about 2.) 
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Fic. 3. Two Morse curves for HF. The F atom centers are 
2.25A apart as in FHF-. D,=6.4 ev. 


1.15A which is, however, the distance from the 
center of the ion to a fluorine nucleus, because 
it is assumed in the application of the theory of 
small vibrations that the hydrogen atom resides 
midway between the Suorine nuclei and nothing 
in this calculation indicates two potential 
minima. Hence the half-distance between the 
fluorine atoms is 1.13A from crystal structure 
study and 1.15A from spectroscopy via the 
theory of small vibrations. The average of i.14A 
will be taken for the next consideration. Since 
the distance from the center of the ion to the 
potential minimum has been estimated to be 
0.26A, the H—F distance is 0.88A. This value 
does not seem unreasonable when compared with 
the H—F distance in monomeric hydrogen 


of 


are 


fluoride, 0.92A, obtained from its band spectro- 
scopy.® Another way of making the comparison 
is to say that the distance (/) from the ion center 
to the potential minimum would be 0.22A 
(1.14—0.92) if the proton remains at the same 
equilibrium distance from the fluorine nucleus 
in the FHF~ ion as in HF. These various pos- 
sibilities are summarized in Table I. It appears 
that these values check sufficiently well to 
warrant the belief that the doublet under con- 
sideration is due to the ‘‘tunnel effect.” It is to 
be expected in FHF~- and a further calculation 
of the potential hill also favors this view. The 
height of the potential barrier can be obtained 
from 
= = 38,750 cal./mole, 


with »vp=3719.7 cm and q=2.7. The zero- 
point energy is one-half the fundamental fre- 
quency v3 multiplied by Planck’s constant and 
amounts to 5316 calories. The barrier has a value 
of 33,400 calories. This quantity can be com- 
pared with the potential hill formed when two 
Morse curves for HF are approached in such a 


°G. Herzberg, Molecular Spectra and Molecular Structure 
(Prentice-Hall, Inc., New York, 1939), p. 488. 
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TABLE I. 


Distance Remarks 


a 3(F-F) 1.13A Crystal structure of 

b H-F 1.15 k=0.4 megadyne/cm from 
3719.7 cm and Badger’s rule 
(Really (F—F) distance!) 

Average: a and b 


c 43(F-F) 1.14 

d H-F 0.92 In HF monomer (band spectra) 
e l 0.22 c¢ minus d 

f l 0.26 Tunnel effect 

g H-F 0.88 HF distance in FHF-; ¢ minus f. 


way as to make the F—F distance 2.25A (see 
Fig. 3). The height calculated is 38,065 calories/ 
mole. For the Morse curve!® the data on HF 
given by Herzberg® were used. The heat of 
dissociation of HF is 6.4 ev and w,.=4141.305 
cm~!, The zero-point energy is calculated to be 
5918 cal./mole. Hence the potential hill inter- 
vening is 32,150 calories/mole. This height is of 
the same order as the barrier in FHF~ as was 
to be expected. The problem of the FHF™ ion 
has been discussed by Pauling." 

10 Reference 9, pp. 108-9. 

uL,. Pauling, The Nature of the Chemical Bond (Cornell 


University Press, Ithaca, New York, 1940), second edition, 
chapter 9. 
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An investigation was made of the errors resulting from electrode polarization in measurements 
of the dielectric constants of solutions containing electrolytes. The applicability of two methods, 


one used by Fricke and Curtis and the other by Oncley, to eliminate these errors, was con- 
sidered. There is introduced a modification of Oncley’s method, which was found applicable at 
higher electrolyte concentrations than Oncley’s method. The relative effectiveness of the 
modified method and that of Fricke and Curtis was tested by employing the two methods to 
correct experimental data obtained for aqueous solutions of glycine containing various amounts 
of electrolytes. The measurements were made at frequencies ranging from 10 to 5000 kc, and at 
conductivities ranging from 0.29 10-4 to 1.7410-* cm™. The corrected values of 
the dielectric constant were compared with known values in the literature and found to be 
in good agreement. The limitations of the methods are considered and the usefulness of the 


I. INTRODUCTION 
HE polarization of dielectric constant cell 
electrodes is an important source of error 
in the determination of the dielectric properties 


methods in cases where the dielectric constant varies with frequency is discussed. 


of dilute electrolyte solutions. This phenomenon 
becomes increasingly important as the frequency 
is lowered and as the conductivity of the solution 
is raised. Where precise dielectric constant 
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610 THOMAS 
determinations are desired steps must be taken 
either to reduce the magnitude of the polarization 
to the point where it may be neglected or 
appropriate correction measures must be em- 
ployed. It is a common practice to reduce the 
amount of polarization by using platinized 
electrodes and for many purposes this procedure 
constitutes an acceptable solution to the diffi- 
culty. 

For measurements on biological systems 
platinized electrodes have several disadvantages. 
They require frequent rechecking because of the 
instability of the platinization and they may 
effect, catalytically, a change in the system under 
investigation.! Furthermore, even though plati- 
nized electrodes are used, for precise measure- 
ments at low frequencies electrode polarization 
may still be significant. For these reasons the 
elimination of electrode polarization errors for 
bare platinum electrodes was investigated. 

Aqueous solutions of glycine, containing 
known amounts of electrolytes, were studied. 
Two methods were employed to correct the 
experimental data for electrode polarization. 
As a test of the effectiveness of the polarization 
correction methods the dielectric constants 
calculated from the experimental measurements 
were compared with known values. 


II. ELECTRODE POLARIZATION CORRECTION 
METHODS 

The problem of correcting dielectric cell 
capacitance measurements for electrode polari- 
zation capacity is complicated by the fact that 
the frequency dependence of the polarization 
capacitance is not definitely known. Although 
it varies with the kind of electrodes and with 
the solution used, the polarization capacitance 
C,, for limited ranges of frequency, can usually 
be expressed by an equation of the form 


Cp,>=Mr™, (1) 


where vy is the frequency, and M and m are 
empirical constants.? 

When the capacitance of a dielectric cell is 
measured it is usually represented electrically in 
terms of an equivalent parallel capacitance C. 
This capacitance is thus a measure of both the 


1J. W. Mehl and C. L. A. Schmidt, J. Gen. Physiol. 18, 
467 (1935). 
21. Wolfe, Phys. Rev. 27, 755 (1926). 
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electrode polarization capacitance and the ca- 
pacitance of the solution contained between the 
cell electrodes. In general the dielectric constant 
of the solution depends on the frequency of 
measurement,? and this of course causes a 
corresponding dependence of the capacitance of 
the solution contained between the cell elec- 
trodes. Since the polarization capacitance of the 
electrodes is also frequency dependent, as Eq. (1) 
indicates, the variation of the cell capacitance 
with changes in frequency may be exceedingly 
complex. | 

Two methods have been used to treat the 
problem of electrode polarization effects in the 
measurement of dielectric constants of solutions. 
The first method to be discussed was described 
by Fricke and Curtis,‘ and was used in a study 
of the dielectric properties of suspensions. In 
their work platinized electrodes were used but 
the method, in principle, is applicable to either 
bare or platinized electrodes. The second method 
was described by Oncley® and was used in studies 
of the dielectric properties of protein solutions 
using a cell with sandblasted platinum electrodes. 

In the method described by Fricke and Curtis 
measurements of the capacitance of the dielectric 
cell are made for two different spacings of the 
cell electrodes. When this is done the dielectric 
constant of the solution can be obtained without 
knowing the absolute value of the polarization 
capacitance. The latter appears only as a 
parameter in the calculations, and is eliminated. 

The dielectric constant ¢ of the solution is 
calculated from the measured capacitances C; 
and C2, corresponding to the electrode spacings 
d, and dz, by means of the equation® 


e=[d2?(C2— Co) —dy?(Ci— Co) J/LK(d2—di) J. (2) 


Here Cy and K are constants for the particular 
dielectric cell. Their significance and evaluation 
are discussed in Section III of this paper. 
Oncley’s method, referred to previously, makes 
use of an empirical equation to correct for 
electrode polarization. The dielectric constant is 


3P. Debye, Polar Molecules (Chemical Catalogue Com- 
pany, New York, New York, 1929). 
a > i. Fricke and H. J. Curtis, J. Phys. Chem. 41, 729 

5 J. L. Oncley, J. Am. Chem. Soc. 60, 1115 (1938). 

6 A derivation of this equation and a statement of the 
conditions under which it is valid are given in the Appendix 
to this paper. 
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Fic. 1. Photograph of micrometer type cell showing lead 
and switching arrangement. 


calculated by means of the equation,® 
(3) 


where C is the observed cell capacitance for an 
electrode spacing d, Cy and K are constants of 
the cell, and AG?*y~ is the correction term for 
polarization. In the correction term, G is the 
cell conductance, v is the frequency, and A isa 
constant for the solution under consideration. 
The value of A is the slope of the straight line 
obtained when (C—Cp) is plotted against 
in a frequency region where the dielectric 
constant does not vary. When measurements are 
extended to a region where the dielectric constant 
varies with frequency, it is assumed that the 
same value of A still prevails. 

Oncley found that the polarization term 
AG*r> allowed accurate corrections for electrode 
polarization to be made only for frequencies 
greater than 25 kc. For lower frequencies it was 
found that the simple correction term AG?y— was 
not valid and the curves were no longer linear. 

In a recent paper Smiley and Smith’ have 


7 W. G. Smiley and A. K. Smith, J. Am. Chem. Soc. 64, 
624 (1942), 


reported a study of electrode polarization in the 
frequency range 1 to 30 kc. These authors 
conclude from their studies, as reported by 
Ferry and Oncley,® that at frequencies less than 
25 kc the correction term AG*v~} is not applicable. 

To meet this difficulty at low frequencies, 
Ferry and Oncley® assumed that the polarization 
correction could be considered as the sum of 
two terms: one, a function of the conductivity 
of the solution, and the other, a function of the 
nature of the solute. The difference in polariza- 
tion of two solutions of equal conductivity was 
then found proportional to v}. Accordingly, 
then, polarization corrections could be made by 
comparing the unknown solution with a known 
salt solution of equal conductivity. With this 
method, Oncley and Ferry have made corrections 
for polarization in protein solutions with con- 
ductivities ranging up to 4010-* cm 
at frequencies as low as 10 kc. 

It will be shown later, in the present paper, 
that Oncley’s original correction scheme can be 
retained and the use of a comparison solution 
avoided, even at frequencies as low as 10 kc and 
for conductivities of the order 200X10-* ohm~ 
cm, 

Ill. EXPERIMENTAL 
Apparatus 


The cell used in this investigation is shown in 
Fig. 1. This cell is of the micrometer type with 
two circular platinum sheet electrodes, one 
mounted on the spindle of a micrometer screw, 
and the other fixed. This type of cell is especially 
useful for making electrode polarization correc- 
tions according to Fricke’s method. The electrode 
spacing may be varied from 0.05 to 10.0 mm 
with a precision of setting equal to +0.002 mm. 
The temperature of the cell was controlled to 
within +0.1°C by circulating mineral oil from 
a thermostat through the glass jacket enclosing 
the cell. The volume of solution ordinarily used 
is 2 ml. 

By mounting the cell on the condenser used as 
a capacitance standard, as shown in Fig. 1, 
several advantages were derived. Short direct 
leads resulted in a low value of lead inductance, 
and a rigid and easily reproduced cell mounting 
resulted in constant lead capacitances. 


( oly Ferry and J. L. Oncley, J. Am. Chem. Soc. 63, 272 
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Fic. 2. Schematic diagram of bridge circuit arranged for 
substitution measurements. Components used in set-up 
are: Bridge, General Radio Type 516-CS1; Oscillator, 
General Radio Type 684A; Detector, Sargent Radio Model 
11VA and Dumont Oscilloscope Type 168 equipped with 
General Radio Type 814 Amplifier and 8 Type 814P 
filters ; Capacitance Standard, General Radio Precision Con- 
denser, Type 722D; Resistance Standard, General Radio 
Decade Resistance Box equipped with inductance compen- 
sating network, Rc, Cc (see reference 8). Cell, see Fig. 1. 


A bridge substitution method was used to 
determine the capacitance and the resistance of 
the cell. At each of the various frequencies 
employed, a precision variable condenser was 
substituted for the cell by means of the switch S, 
indicated in Fig. 2, which is a schematic diagram 
of the circuit used. After the substitution the 
bridge was rebalanced by adjustment of the 
precision condenser and an inductance-compen- 
sated resistance box. The various items of 
equipment making up the circuit have been 
identified on the diagram. It should be noted 
that in the electrical set-up employed in this 
paper all lead inductances are measured from 
the terminals of the precision condenser. The 
latter may then be considered as the reference 
terminals of the system. 

The use of a substitution method effectively 
eliminates from the measurements errors which 
arise from the bridge proper. The further errors 
which must be considered are those due to the 
inductances L, Ls, Lr, of the leads between the 
capacitance and resistance standards and the 
cell. These inductances affect both the cell 
resistance and capacitance obtained by the 
substitution method. In the work reported in 
this paper only the low frequency measurements 


9J. L. Oncley, J. Phys. Chem. 44, 1103 (1940). 
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of resistance are used, for which lead inductance 
effects are negligible. The following discussion is 
therefore directed to the elimination of errors in 
the measured capacitances. 

The effect of lead inductances on a capacitance 
determination may be represented by the 
equation 


(4) 


where Cs is the capacitance setting of the 
precision condenser, C is the true terminal 
capacitance of the cell, and the delta terms are 
capacitance errors to be described below. 

The series inductance L of the cell contributes 
capacitance errors, AC, and AC;. This inductance 
is in series with the unknown capacitance C of 
the cell and gives rise to the capacitance incre- 
ment AC,=[C—C/(1—w*LC)] for a frequency 
v=w/2n. Further, since L is in series with the 
resistance R of the cell, it gives rise to a capaci- 
tance increment AC,=(—L/R?). The equations 
expressing the delta terms are valid when 
R°>.*L, which condition holds for the fre- 
quencies and resistance values involved here. 
A detailed treatment of lead errors is given by 
Davies.” 

The inductance of the standard condenser Ls 
behaves in a manner similar to L. At high 
frequencies the capacitance of the standard 
condenser is increased by an amount AC, 
=[Cs—Cs/(1—w*L sCs) ]. Because the standard 
condenser has a high parallel resistance, the 
correction AC,, due to the condenser resistance 
being in series with Ls, is very small and can be 
completely neglected in the work discussed here. 

The inductance Lr of the standard resistor 
Rs gives rise to an effective capacitance 
C.=(—Lr/Rs?*). In making a substitution meas- 
urement of the cell capacitance, the initial 
bridge balance is obtained with the standard 
resistance set at its maximum value, equal to 
260.2 ohms. After the cell is removed from the 
circuit and the standard condenser substituted 
for it, Rs is reduced from its maximum value to 
another value, which depends on the cell 
resistance. Hence the effective capacitance C, of 
the standard resistor changes during a given cell 
capacitance measurement and may cause large 


10R. M. Davies, Phil. Mag. 20, 75 (1935). 
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errors if ignored. This effective capacitance 
change AC. was reduced to less than +0.7 ppf 
for frequencies between 5000 and 10 kc by 
means of a modification of the inductance 
compensation scheme described by Oncley.® A 
description of the method used to evaluate the 
capacitance errors in the resistance standard is 
being prepared for publication. 

For the cell shown in Fig. 1, Z was found to be 
0.1110-® henry. This value was calculated 
from measurements of the variation in effective 
capacitance of the cell at high frequencies. This 
value of L is sufficiently small that the capaci- 
tance AC, may be neglected for all values of the 
cell resistance R used in the present measure- 
ments. For a minimum value of R equal to 450 
ohms, the maximum value of AC, is approxi- 
mately 0.5 wuf. This value is slightly less than 
the maximum capacitance error of the standard 
resistor Rs, mentioned above, and it has there- 
fore been neglected. 

The capacitances AC, and AC, remain to be 
considered. These were caused to cancel each 
other effectively by adjusting L, to equal L.° 


Calibration of the Cell 


Calibration measurements were made with 
water and dilute aqueous potassium chloride 
solutions as standards of known dielectric 
constant. A frequency of 1000 kc was used in 
order to avoid electrode polarization. Measure- 
ments were made of the parallel capacitance of 
the cell for a series of electrode spacings between 
0.0200 and 0.6000 cm. The cell capacitance C is 
given by the equation" 


C=+Ke/d, (5) 


where Cy is the distributed capacitance, K is 
approximately equal to the electrode area divided 
by 4, ¢« is the dielectric constant of the cell 
contents, and d is the electrode spacing. The 
constants Cy and K were evaluated from straight 
line graphs of C plotted as function of 1/d. The 
average values obtained were Cyo=3.5+0.5 upf 
and K=0.1700+0.0005 yuf cm. They were used 
in reducing the capacitance measurements, for 
the glycine solutions, to dielectric constants as 
discussed subsequently. 


1K. S. Cole and H. J. Curtis, Rev. Sci. Inst. 8, 333 
(1937). 
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Materials 


Glycine was chosen for the dielectric constant 
measurements for two reasons. First, the di- 
electric constants of aqueous glycine solutions 
are known with considerable accuracy” for all 
obtainable concentrations of the material. They 
serve, therefore, as convenient standards to 
evaluate the adequacy of the methods for making 
polarization corrections. Second, glycine solu- 
tions were to be used as solvents in a study of 
the dielectric properties of lactoglobulin in this 
laboratory over the same frequency range 
employed in the present investigation and it was 
desirable, therefore, to know the behavior of the 
cell when filled with glycine solutions. 

The glycine was obtained from the Eastman 
Kodak Company and was recrystallized once 
from 90 percent methanol. The solutions studied 
contained 0.5, 1.5, and 2.5 moles of glycine per 
liter. In order to vary the conductivity of the 
solutions, they were made up to contain from 
0.0001 to 0.0015 mole per liter of potassium 
chloride. In this way the conductivity of the 
solutions was caused to vary from 29X10-® to 
10-* ohm™ 


IV. RESULTS 
Tables I, II, and III summarize the results 


for nine solutions of glycine. For the sake of 


TaBLe I, Summary of cell capacitances and dielectric 
constants for glycine solutions. Electrode polarization cor- 
rections made by the Fricke and Curtis method. 


1 2 3 |4|5 6 7 8 9 10 
Observed capacitance, 
Dielectric 
| constant 
Electrode spacing, d € 
d= 
d =0.123 cm | 0.173 d =0.273 cm (Eq. (2)] 
cm 
Fre- 
quency} Solution® 
inke | D E F D E F D E F 
5000 | 158.4 158.5|112.6|71.4 71.7 71.8|114.7 115.5 116.8 
1000 | 158.2 159.0}113.5|71.1 71.5 71.3|113.9 114.6 113.0 
500 | 158.8 158.3} 113.6| 71.1 71.5 71.9}113.6 115.0 115.5 
100 | 159.5 160.7 | 115.9} 71.2 72.0 72.7/|113.5 115.1 114.9 
160.6 73.2 74.6/113.6 114.9 113.8 
40 | 161.5 171.2 |124.8|]71.7 74.1 75.9|113.7 115.0 113.3 
30 | 163.2 179.6] 131.4|72.0 75.7 78.6|113.6 114.7 113.5 
20 | 166.8 201.7|149.0|}72.9 804 85.9!/114.1 115.3 114.5 
10 | 185.7 306.3 | 232.0} 76.7 102.3 120.2|114.0 117.2 118.8 


*Solutions D, E, F contain 1.5 M glycine and have conductivities 
of 55.7, 139, and 174 X10-§ ohm~! cm“, respectively. Conductivities 
were brought to these values by the addition of KCl. 


12 J. Wyman, Jr. and T. L. McMeekin, J. Am. Chem. Soc. 
55, 915 (1933). 
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BRIDGE 
CELL) 
R RM 
CAPACITANCE 
STANDARD 
RESISTANCE _! 
STANDARD 


Fic. 2. Schematic diagram of bridge circuit arranged for 
substitution measurements. Components used in set-up 
are: Bridge, General Radio Type 516-CS1; Oscillator, 
General Radio Type 684A; Detector, Sargent Radio Model 
11VA and Dumont Oscilloscope Type 168 equipped with 
General Radio Type 814 Amplifier and 8 Type 814P 
filters ; Capacitance Standard, General Radio Precision Con- 
denser, Type 722D; Resistance Standard, General Radio 
Decade Resistance Box equipped with inductance compen- 
sating network, Rc, Cc (see reference 8). Cell, see Fig. 1. 


A bridge substitution method was used to 
determine the capacitance and the resistance of 
the cell. At each of the various frequencies 
employed, a precision variable condenser was 
substituted for the cell by means of the switch S, 
indicated in Fig. 2, which is a schematic diagram 
of the circuit used. After the substitution the 
bridge was rebalanced by adjustment of the 
precision condenser and an inductance-compen- 
sated resistance box.’ The various items of 
equipment making up the circuit have been 
identified on the diagram. It should be noted 
that in the electrical set-up employed in this 
paper all lead inductances are measured from 
the terminals of the precision condenser. The 
latter may then be considered as the reference 
terminals of the system. 

The use of a substitution method effectively 
eliminates from the measurements errors which 
arise from the bridge proper. The further errors 
which must be considered are those due to the 
inductances L, Ls, Lr, of the leads between the 
capacitance and resistance standards and the 
cell. These inductances affect both the cell 
resistance and capacitance obtained by the 
substitution method. In the work reported in 
this paper only the low frequency measurements 


9J. L. Oncley, J. Phys. Chem. 44, 1103 (1940). 


of resistance are used, for which lead inductance 
effects are negligible. The following discussion is 
therefore directed to the elimination of errors in 
the measured capacitances. 

The effect of lead inductances on a capacitance 
determination may be represented by the 
equation 


(4) 


where Cs is the capacitance setting of the 
precision condenser, C is the true terminal 
capacitance of the cell, and the delta terms are 
capacitance errors to be described below. 

The series inductance L of the cell contributes 
capacitance errors, AC, and AC. This inductance 
is in series with the unknown capacitance C of 
the cell and gives rise to the capacitance incre- 
ment AC,=[C—C/(1—w*LC)] for a frequency 
v=w/2n. Further, since L is in series with the 
resistance R of the cell, it gives rise to a capaci- 
tance increment AC,=(—L/R?). The equations 
expressing the delta terms are valid when 
R’*>w*L, which condition holds for the fre- 
quencies and resistance values involved here. 
A detailed treatment of lead errors is given by 
Davies.!° 

The inductance of the standard condenser Ls 
behaves in a manner similar to L. At high 
frequencies the capacitance of the standard 
condenser is increased by an amount AC, 
=[Cs—Cs/(1—w*LsCs) ]. Because the standard 
condenser has a high parallel resistance, the 
correction AC,, due to the condenser resistance 
being in series with Ls, is very small and can be 
completely neglected in the work discussed here. 

The inductance Le of the standard resistor 
Rs gives rise to an effective capacitance 
C.=(—Lr/Rs*). In making a substitution meas- 
urement of the cell capacitance, the initial 
bridge balance is obtained with the standard 
resistance set at its maximum value, equal to 
260.2 ohms. After the cell is removed from the 
circuit and the standard condenser substituted 
for it, Rs is reduced from its maximum value to 
another value, which depends on the cell 
resistance. Hence the effective capacitance C, of 
the standard resistor changes during a given cell 
capacitance measurement and may cause large 


10R. M. Davies, Phil. Mag. 20, 75 (1935). 
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errors if ignored. This effective capacitance 
change AC. was reduced to less than +0.7 pyf 
for frequencies between 5000 and 10 kc by 
means of a modification of the inductance 
compensation scheme described by Oncley.® A 
description of the method used to evaluate the 
capacitance errors in the resistance standard is 
being prepared for publication. 

For the cell shown in Fig. 1, Z was found to be 
0.11X10-* henry. This value was calculated 
from measurements of the variation in effective 
capacitance of the cell at high frequencies. This 
value of L is sufficiently small that the capaci- 
tance AC, may be neglected for all values of the 
cell resistance R used in the present measure- 
ments. For a minimum value of R equal to 450 
ohms, the maximum value of AC, is approxi- 
mately 0.5 wuf. This value is slightly less than 
the maximum capacitance error of the standard 
resistor Rs, mentioned above, and it has there- 
fore been neglected. 

The capacitances AC, and AC, remain to be 
considered. These were caused to cancel each 
other effectively by adjusting L, to equal L.® 


Calibration of the Cell 


Calibration measurements were made with 
water and dilute aqueous potassium chloride 
solutions as standards of known dielectric 
constant. A frequency of 1000 kc was used in 
order to avoid electrode polarization. Measure- 
ments were made of the parallel capacitance of 
the cell for a series of electrode spacings between 
0.0200 and 0.6000 cm. The cell capacitance C is 
given by the equation" 


C=Co+Ke/d, (5) 


where Co is the distributed capacitance, K is 
approximately equal to the electrode area divided 
by 42, € is the dielectric constant of the cell 
contents, and d is the electrode spacing. The 
constants Cy and K were evaluated from straight 
line graphs of C plotted as function of 1/d. The 
average values obtained were Cyo=3.5+0.5 upf 
and K=0.1700+0.0005 puuf cm. They were used 
in reducing the capacitance measurements, for 
the glycine solutions, to dielectric constants as 
discussed subsequently. 


1K. S. Cole and H. J. Curtis, Rev. Sci. Inst. 8, 333 
(1937), 
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Materials 


Glycine was chosen for the dielectric constant 
measurements for two reasons. First, the di- 
electric constants of aqueous glycine solutions 
are known with considerable accuracy” for all 
obtainable concentrations of the material. They 
serve, therefore, as convenient standards to 
evaluate the adequacy of the methods for making 
polarization corrections. Second, glycine solu- 
tions were to be used as solvents in a study of 
the dielectric properties of lactoglobulin in this 
laboratory over the same frequency range 
employed in the present investigation and it was 
desirable, therefore, to know the behavior of the 
cell when filled with glycine solutions. 

The glycine was oLtained from the Eastman 
Kodak Company and was recrystallized once 
from 90 percent methanol. The solutions studied 
contained 0.5, 1.5, and 2.5 moles of glycine per 
liter. In order to vary the conductivity of the 
solutions, they were made up to contain from 
0.0001 to 0.0015 mole per liter of potassium 
chloride. In this way the conductivity of the 
solutions was caused to vary from 29X10-* to 
174X10-* cm—. 


IV. RESULTS 
Tables I, II, and III summarize the results 


for nine solutions of glycine. For the sake of 


TaBLe I, Summary of cell capacitances and dielectric 
constants for glycine solutions. Electrode polarization cor- 
rections made by the Fricke and Curtis method. 


1 2 3 |4|5 6 7 8 9 10 
Observed capacitance, puf 
on Dielectric 
(C —Co) | (C —Co) 
Electrode spacing, d € 
d= 
d =0.123 cm | 0.173 d =0.273 cm {Eq. (2)] 
cm 
Fre- 
quency, Solution® 
in ke D E F D E F D E F 
5000 | 158.4 158.5|112.6|71.4 71.7 71.8|114.7 115.5 116.8 
1000 | 158.2 159.0/113.5|71.1 71.5 71.3}113.9 114.6 113.0 
500 | 158.8 158.3|113.6|71.1 71.5 71.9/113.6 115.0 115.5 
100 | 159.5 160.7|115.9| 71.2 72.0 72.7|113.5 115.1 114.9 
50 | 160.6 166.9/121.3|71.5 73.2 74.6|113.6 114.9 113.8 
40 | 161.5 171.2 |124.8|71.7 74.1 75.9|113.7 115.0 113.3 
30 | 163.2 179.6|131.4|72.0 75.7 78.6|113.6 114.7 113.5 
20 | 166.8 201.7|}149.0}72.9 80.4 85.9/114.1 115.3 114.5 
10 | 185.7 306.3 | 232.0| 76.7 102.3 120.2|114.0 117.2 118.8 


*Solutions D, E, F contain 1.5 M glycine and have conductivities 
of 55.7, 139, and 174 X10-* ohm~! cm™, respectively. Conductivities 
were brought to these values by the addition of KCI. 


12 J. Wyman, Jr. and T. L. McMeekin, J. Am. Chem. Soc. 
55, 915 (1933). 
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TABLE II. Summary of corrected capacitances and dielectric constants for glycine solutions. Polarization corrections 
according to modification of the Oncley method. 


1 2 3 | 8 6 7 8 9 10 | 11 12 13 
Commies Dielectric constant, ¢ 


Electrode spacing, d 


d =0.123 cm d=0.173 cm 


Electrode spacing, d 


d=0.123 cm d=0.173 cm 


D E E 


E E F 


158.4 158.5 113.3 114.6 114.7 115.3 114.6 

1000 158.1 159.0 113.1 113.8 113.4 114.4 115.0 115.1 115.8 115.4 
500 158.7 158.3 113.0 113.8 113.4 114.8 114.5 115.0 115.8 115.4 
100 159.0 157.8 113.2 113.7 113.5 115.0 114.2 115.2 115.7 115.5 
50 158.8 157.6 113.3 113.4 113.6 114.9 114.0 115.3 115.4 115.6 
40 158.9 157.7 113.3 113.4 113.6 115.0 114.1 115.3 115.4 115.6 
30 159.0 157.5 113.3 113.2 113.1 115.0 114.0 115.3 115.2 115.1 
20 159.0 157.8 112.7 113.6 115.0 114.2 115.5 114.7 115.6 


158.7 163.3° 113.9 


113.5 


114.8 118.1° 115.5 116.0 115.5 


brought to these values by the addition of KCl. 
b Measurements not made for solution F at this electrode s 


brevity, detailed data are presented only for 
the solutions of the 1.5 M glycine. 

In Table I, columns 2-7, are shown observed 
capacitances for three solutions of 1.5 M 
glycine of different conductivities and for two 
values of electrode spacing. In columns 8-10 
appear the dielectric constants of the solutions 
calculated from the values of C by means of 
Eq. (2) according to the method of Fricke and 
Curtis. The capacitance values in columns 2-7 
show clearly the regular increase due to electrode 
polarization that occurs with decreasing fre- 
quency and increasing conductivity. For the 
solutions of highest conductivity the polarization 
becomes so large that at low frequencies the 
observed capacitances are approximately twice 
the high frequency values. At the three highest 
frequencies no polarization effects are apparent. 
The capacitance variations that occur at these 
frequencies are within the limits of experimental 
error. 

The dielectric constants shown in columns 
8-10 should be free from electrode polarization 
effects if the Fricke and Curtis method is valid. 
The agreement manifest between the values of 
¢ for solutions of different conductivity at high 
and low frequencies shows this to be the case. 
No significant trends with frequency or con- 
ductivity are apparent, the values of ¢« for the 
lower frequencies and higher conductivities being 
in good agreement with the values for the higher 


° This value is probably in error; not included in averages = padng III. 


®Solutions D, E, F contain 1.5 M glycine and have conductivities of 55.7, 139, and 174 X10-* ohm™~! cm™|, respectively. Conductivities were 


TABLE III. Comparison of dielectric constants of glycine 
solutions showing effectiveness of methods used to correct 
for electrode polarization. 


1 | 2 3 oy @ 
Solution Dielectric constant, « 
— Results of present work 
tivity Fricke and Curtis | Modified Oncley 
at 10 kc method method 
Glycine | |Accepted Average Average 
moles/liter | 10° | valued | Mean* | devia- | Means | devia- 
tion tion 
92.2> | +0.06 
(A) 0.5 29.0 90.9 +0.11 91.4 +0.07 
91.7 +0.16 
(B) 0.5 75.4 91.74 91.5 +0.27 92.1 +0.14 
91.6 +0.43 
(C) 0.5 135 91.9 +0.19 91.8 +0.30 
114.8 +0.17 
(D) 1.5 55.7 .| 113.8 +0.28 115.2 +0.18 
114.3 +0.26 
(E) 1.5 139 114.52 115.2 +0.48 115.5 +0.31 
— 6 — 
(F) 1.5 174 114.9 +1.42 115.3 +0.27 
136.7 +0.08 
(G) 2.5 60.4 139.3 +0.61 137.4 +0.11 
138.1 +£0.06 
(H) 2.5 99.7 137.30 138.2 +0.24 138.3 +0.24 
137.2 +0.44 
q@) 2s 140 138.1 +1.07 138.5 +0.22 


® Mean of 9 values over frequency range 10 to 5000 kc. 
b In this column the upper figure is for d =0.123 cm; lower figure is 
for d =0.173 cm. 
¢ No measurements were made for d =0.123 cm. 
4 Values for dielectric — from J. Wyman, Jr., Phys. Rev. 35, 
623 (1930) PS gp J. an, Jr., and T. L. McMeekin, J. Am. 
Chem. Soc. 55, 913 


frequencies and lower conductivities, where 
electrode polarization is negligible. Average 
values from columns 8-10 of Table I are sum- 
marized in Table III, together with similar 
data for 0.5 and 2.5 M glycine, for comparison 
with the results obtained by an application of a 


| 
Solution® 
Frequency 
in ke i D 


modification of Oncley’s method to the same 
experimental data. 

In Table II are presented the results obtained 
when a modification of Oncley’s method is used 
to correct the experimental data for electrode 
polarization. Oncley’s method was modified in 
order to be able to make corrections at low 
frequencies without introducing the additional 
complication of the comparison method that 
Ferry and Oncley® found necessary. When values 
of the observed capacitance C are plotted as a 
function of G?y-> or simply v5, because G? 
changes very little in the frequency region 
concerned," curves of the type shown in Fig. 3 
are obtained. These curves are linear at all 
except the highest conductivities for frequencies 
down to about 25 kc. Below this frequency the 
curves are no longer linear. This behavior was 
noted by Ferry and Oncley and resulted in 
their adoption of a comparison method, men- 
tioned previously. 

In an effort to avoid the additional complica- 
tion introduced by the comparison method, tests 


8 


CAPACITANCE (MMF) 


2 4 6 
ys x 10” 
Fic. 3. Curves showing relation between cell capacitance 
Cc and frequency v to —1.5 power for 1.5 molar glycine solu- 
tions of varying conductivity o and for different dielectric 
cell electrode spacings d. 
Curve 1, 55.7X10-¢ d=0.123 cm. 
Curve 2, 55.7X107-* d=0.173 cm. 
Curve 3, 91.4X10-* cm, d=0.123 cm. 
Curve 4, c= 91.4 10-6 ohm™ d=0.173 cm. 
Curve 5, ¢=139 d=0.123 cm. 
Curve 6, ¢=139 ohm™ d=0.173 cm. 
Curve 7, ¢=174 ohm™ d=0.173 cm. 


‘3 In all the experiments here reported the total change in 
conductance amounts to a decrease of 1 percent or less for 
the frequency range 10 kc to 500 ke. 
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CE (MMF) 
wn 


CAPACITAN! 
8 


1005 4 
yt? x 107 

Fic. 4. Curves showing relation between dielectric cell 

capacitance C and frequency v to —1.7 power for 1.5 molar 

glycine solutions of varying conductivity o and for different 


dielectric cell electrode spacings d. 


Curve 1, d=0.123 cm. 
Curve 2, «= 55.7X10-§ ohm d=0.173 cm. 
Curve 3, c= 91.4X10-§ ohm™ cm™, d=0.123 cm. 
Curve 4, c= 91.4X10-§ d=0.173 cm. 
Curve 5, ¢=139 d=0.123 cm. 
Curve 6, c=139 X10-§ d=0.173 cm. 
Curve 7, ¢=174 X10-§ d=0.173 cm. 


were made in the present work to see if the 
adoption of a different frequency dependence 
would yield straight lines. By successive trials 
it was found that when values of C were plotted 
as a function of v—’, straight lines resulted. 
This behavior was found not only for the glycine 
solutions, for which data are reported here, but 
in addition for aqueous solutions of potassium 
chloride with concentrations ranging from zero 
to 0.002 N. Typical curves of C versus v™? are 
shown in Fig. 4. 

The fact that graphs of C versus v—-’ are linear 
would appear to indicate that the constant m in 
Eq. (1) is equal to 0.3" for the lower frequencies, 
at least, rather than 0.5 as was assumed by 
Oncley.> That m varies with the frequency is 
well known. For example, Wolfe? reported m to 
be approximately 0.5 in the upper portion and 
approximately 0.3 in the lower portion of a 
frequency range of 0.2 kc to 200 kc. The present 
results, however, seem to indicate that m is not 
critical for high frequencies and for low conduc- 
tivities, and that m=either 0.5 or 0.3 represents 
the data adequately. At the very low frequencies, 


14 See Appendix 2. 
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however, as a comparison of Figs. 3 and 4 
shows, the smaller value of m is essential. 

The values of C’, the corrected capacitance of 
the cell, which appear in columns 2-7 of Table II, 
were therefore obtained according to the equation 


C’=(C—Cy— Sv"). (6) 


The constant S was obtained from the straight 
lines resulting when C was plotted against y™’. 
The values of C’ are seen to be independent of 
frequency and conductivity. The excellent agree- 
ment between high and low frequency values 
confirms the use of Eq. (6) for the elimination 
of the polarization effects which occur at the 
lower frequencies. 

The dielectric constants which appear in 
columns 8-13 of Table II were calculated 
according to the relation e=C’d/K. 

Averages of the values in columns 8-13 have 
been assembled in Table III for comparison 
with the results obtained by the method of 
Fricke and Curtis. 

An inspection of Table III reveals that the 
Fricke and Curtis method and the modified 
Oncley method are equally adequate to eliminate 
polarization effects. The dielectric constants 
obtained by an application of these methods to 
the experimental data are seen to agree well with 
the values of Wyman and McMeekin.” This is 
true for solutions with conductivities ranging 
from 0.29 to 1.74X10-* No de- 
pendence of the dielectric constant on the 
conductivity of the solution is apparent. The 
differences between the ‘‘mean”’ and “‘accepted”’ 
values is about the same for the two methods. 
With one exception these differences are of an 
order of magnitude consistent with the uncer- 
tainty in the cell constants. 

The quantities listed as ‘‘average deviations” 
in columns 5 and 7 represent deviations among 
the values for different frequencies. In most 
instances these deviations are much less for the 
modified Oncley method than for the method of 
Fricke and Curtis. They probably reflect errors 
arising from incomplete compensation of the 
lead inductance in the standard resistor, as 
mentioned earlier. The greater deviation in the 
case of the results with the Fricke and Curtis 
method is to be expected since twice as many 
observations, each subject to the same uncer- 


THOMAS M. 


SHAW 


tainty, enter into the calculations, as a com- 
parison of Eqs. (2) and (6) shows. 


V. CONCLUSIONS 


The data presented demonstrate that errors 
due to electrode polarization in dielectric con- 
stant measurements may be eliminated equally 
well by the application of the Fricke and Curtis 
method or the modification of Oncley’s method, 
described in the present paper. This conclusion 
applies, of course, only for the frequency range 
and conductivity range investigated. There 
appears to be no reason why the Fricke and 
Curtis méthod cannot be extended to much 
lower frequencies and to higher conductivities. 
On the other hand the extension of the modified 
Oncley method is questionable. Since the fre- 
quency dependence of the polarization is known 
to vary with the frequency, it is quite likely that 
for a precise correction of the polarization at 
lower frequencies it would be necessary to take 
this fact into account. 

Of the two methods for eliminating polariza- 
tion errors, there is nc doubt of the more general 
utility of the Fricke and Curtis method. This is 
true despite the fact that the dielectric constants 
obtained with their method are less precise than 
the values obtained with the modification of 
Oncley’s method. The greater utility arises from 
the fact that the Fricke and Curtis method can 
be used to eliminate polarization effects in 
measurements made at a single frequency only. 
Oncley’s method or the present modification of 
it, on the other hand, requires measurements for 
at least two frequencies; moreover, these must 
lie outside of the dispersion region of the material 
under investigation. This latter limitation is 
particularly important in studies of proteins or 
other high molecular weight compounds where 
the dielectric constant may show dispersion over 
an extremely wide range of frequencies. In those 
instances where this limitation does not apply, 
Oncley’s method would appear preferable on 
the basis of increased precision and because it 
may be used for measurements made with a 
simple cell having fixed electrodes. 


APPENDIX 1 
Derivation of Eq. (2) 


It is desired to find an equation expressing the dielectric 
constant of a solution in a dielectric constant cell with 


ELECTRODE POLARIZATION 


Fic. 5. Schematic diagram showing relation between 
measured capacitances and resistances and equivalent 
impedances. C;, Ri, measured parallel capacitance and 
resistance of cell for electrode spacing d:; C2, Re, measured 
parallel capacitance and resistance of cell for electrode 
spacing d2; and Co, distributed capacitance of cell. 


Z»=impedance of electrodes; 

Za=impedance of solution defined by electrode 
spacing 

Zp = increase in impedance of solution due to change in 
electrode spacing from d, to do; 

Z,=total impedance of cell for electrode spacing di; 

Z2=total impedance of cell for electrode spacing do. 


Dashed quantities are equivalent series capacitance and 
resistance components of corresponding impedances. 


polarized electrodes in terms of the measured resistance and 
capacitance of the cell. Referring to Fig. 5, where the 
symbols are defined, it can be seen that the impedance of 
the element of solution, Zz, of length (d2:—d,) is given by, 


Zp=Z2—-Z). (7) 


Making use of the well-known relations between the 
capacitance and resistance of series and parallel equivalent 
circuits, substitution for the impedances in Eq. (7) yields, 


= [w*(C2— Co) +w*(C2— Co)?R? 
— [w?(Ci— Co) RP (8) 
and 
(9) 
Now for the conditions, which hold for all measurements 
reported in this paper where the polarization is significant, 
w?(Ci— Co)*RPKA, w®(Co— Co)*R2K 1, 
and 


(10) 


it follows that, 
Ca =[(C2—Co)R2—(Ci— Co) (11) 


(12) 


and 
Re=(R2—R)). 


Equation (11) may be simplified further. Jones and 
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Christian’ have shown that the equivalent series resist- 
ances of a cell of the type used here may be expressed, 

R,=(d:/Ac)+Rp (13) 
and 

R2=(d2/Ac)+Rp, (14) 
where A is the area of the cell electrodes and @ is the specific 
conductivity of the solution. For the conditions (10), Ri= Ri 
and R.= Re». Thus Eqs. (13) and (14) may be substituted in 
Eq. (11) with the result, 


Ca =([d2*(C2— Co) —d2(Ci— Co) (15) 


Now it follows from the theory for a parallel plate con- 


denser that 
Cp=eK/(d2—d,), (16) 


where K is a constant. The elimination of Cg from Eqs. (15) 
and (16) yields finally, 


e = [d2?(C2— Co) — Co) ]/ KK (d2—d). (2) 


This equation also holds for frequencies where electrode 
polarization is negligible but for which the condition ex- 
pressed in (10) is not satisfied. When the electrode polariza- 
tion is negligible, (Ci—C))=K’'/di, 
K"'(d2—d;), where for a given cell and solution K’ and K”’ 
are constants. If these values are substituted in Eqs. (8) 
and (9), then 

Ca=Ci-C2/(Ci— C2). (17) 


This equation may be used interchangeably with Eq. (2) 
when electrode polarization effects are absent. 


APPENDIX 2 


Derivation of Eq. (3) 

The cell is assumed to be represented by a measured 
parallel capacitance (C—Co) and measured parallel con- 
ductance G=(1/R), as Fig. 5 indicates. These quantities 
are related to the series polarization capacitance Cy of the 
electrodes, and the series capacitance C of the solution in 
the cell by the following equation, 


(C— Cy) = (18) 
which can be obtained by a process similar to that used in 
Appendix 1, for the condition, R*(C— Co)*w<<1. Assuming 
that the polarization capacitance may be represented by 
C,=Mv™ [Eq. 1] by substitution is obtained, 

(C— Co) = (G2v-2-™) M) + (19) 

When R?C%?K1 then C =G?/4x*v°C, which when substi- 
tuted in Eq. (19) yields 

(20) 

where C’ is the cell capacitance corrected for electrode 


polarization. Now, making use of the equation for a parallel 
plate condenser, e= C’d/K, Eq. (20) becomes 


]d/K, (3) 


which, for m=0.5, is the desired equation. 


15 G. Jones and S. M. Christian, J. Am. Chem. Soc. 57, 
272 (1935). 
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Long, heavy lead mirrors are deposited in a reaction tube, and metal is transported there- 
from by free radicals according to the Paneth method. The lead alkyl is collected and the lead 
determined by a micro-method. Assuming that all the alkyl radicals are mirror-active and that 
all the lead is transported as PbR,, the number of radicals can be estimated. The method is 
applied to work on the photolysis of acetone at 2537A. It is shown that the difference in ap- 
parent half-lives of radicals produced at 45°C (1.1107? sec.) and at 100°C (8X 10-3 sec.) may 
be explained on the basis of the instability of free acetyl. The actual concentration of methyl 
radicals in the flowing gas in the experiments reported is estimated to be ~2 X 10-3 mole per- 
cent or 1.7 10-" mole per cc, immediately past the irradiated zone. 


INTRODUCTION 


LTHOUGH methods of absorption spectros- 
copy have been used for determining the 
concentrations of such free radicals as OH! and 
CN,? such methods are not in general available 
for the estimation of the more complicated free 
radicals represented by the alkyls.* Given an 
adequate method of measuring the metal trans- 
ported, the classical mirror methods of Paneth 
and his co-workers‘ would seem to be applicable. 
One method which has been used to count 
alkyl radicals is to deposit a tellurium mirror of 
known weight in a tube and to remove it ac- 
cording to the Paneth method. The time re- 
quired to remove the mirror completely has 
been used to estimate the free radical concentra- 
tion > however, this method is in error for it is 
implicitly assumed in the calculation that every 
radical strikes the mirror at least once in passing 
it, no matter how much the mirror area may be 
reduced. Paneth, Hofeditz, and Wunsch® have 
picked up lead from a large mirror in a similar 
way and have thermally decomposed the lead 
tetramethyl so formed in another portion of 
the reaction tube. After the section containing 
the lead deposit is cut out of the tube, the sample 
may be removed and estimated by turbidimetric 


10. Oldenberg, J. Chem. Phys. 2, 713 (1934). 

2]. U. White, J. Chem. Phys. 8, 459 (1940). Cf. G. 
Kistiakowsky and H. Gershinowitz, ibid., 1, 432 (1933). 

3 Cf. O. Oldenberg, J. Phys. Chem. 41, 293 (1937). 

4See F. O. Rice and K. K. Rice, The Aliphatic Free 
Radicals (Johns Hopkins Press, Baltimore, 1935) for 
references. 

5]. S. A. Forsythe, Trans. Faraday Soc. 37, 312 (1941). 
a 935) Hofeditz and Wunsch, J. Chem. Soc. 372 


methods. Although Paneth and his co-workers 
did not specifically use the method for the esti- 
mation of free radicals, it is interesting to note 
that their paper shows that at least 82 percent 
of the free radicals produced were active in 
removing mirrors; the general inference which 
may be drawn from their work is that very 
probably every methyl radical which strikes a 
mirror reacts with it. In a similar method,’ the 
lead deposited was measured by its radioactivity. 
The latter method is very sensitive, but it is 
not readily susceptible to quantitative tech- 
niques. Furthermore, like the other deposition 
method, it is not quite satisfactory because the 
decomposition of the lead alkyl changes the rest 
of the system under investigation. 

In the modification of the Paneth method 
described in this report, the lead is removed from 
a long heavy mirror by free radicals and col- 
lected as lead alkyl which is then converted to 
lead nitrate and estimated by microanalytical 
techniques. In this way, 25ug, or ~1.210~7 
mole of Pb is conveniently measured. The 
number of free radicals involved may then be 
calculated on the basis of the assumption that 
all the lead is transported as PbR,. The applica- 
tion of this technique is illustrated in work on the 
half-life of radicals from the photolysis of acetone. 


EXPERIMENTS AND DISCUSSION 
Apparatus and Method 
' The apparatus, with the exception of a col- 
lection bulb M, shown in Fig. 1, is in all es- 


7 Burton, Ricci, and Davis, J. Am. Chem. Soc. 62, 265 
(1940). 
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sentials similar to that described in previous 
publications from this laboratory. The same 
careful, high vacuum techniques and methods 
of purifying and handling acetone are employed. 
A departure is in the actual deposition of the 
lead mirrors in the quartz reaction tube: (a) 
the final deposit is much heavier than that in the 
usual mirror technique, and (b) the length of 
the mirror in the tube (it covers the cylindrical 
surface completely and heavily) is fixed to ex- 
ceed a definite amount. In this work, as will be 
shown, 17-20 mm was the necessary minimum 
length. 

The methods of illumination (at 2537A with 
a spiral quartz mercury resonance lamp) and 
of starting and stopping runs were as already 
described.* Temperatures in the irradiated zone 
were controlled manually to within 1.5° with an 
electric fan. Asbestos housing was required to 
give better control at 100°C. The temperature of 
the rest of the reaction tube was between 20° 
and 25°C, depending on the room temperature. 
All experiments were conducted at an input 
pressure of acetone (into the reaction tube) 
of 2.1 mm Hg. Previous calibration showed that 
under such pressure and at room temperature in 
the reaction system, the streaming velocity of 
acetone through the 0.54-cm internal diameter 
quartz tube of length set at ~75 cm was 3.0 
(+.1) 108 cm/sec. 


Collection of Lead 


For the determination of the quantity of 
lead picked up by methyl radicals during the 
course of a run, the condensable products were 
collected in the bulb L kept at liquid nitrogen 
temperature. After isolation from the system 
this mixture was warmed to 35-40°C and dis- 
tilled into a second trap M at —196°C, contain- 
ing about 1 cc of a mixture (in equal volumes) of 
bromine and carbon tetrachloride, previously 
evacuated and stored at this temperature. This 
second trap was then shut off and warmed for 
five minutes at 0-25°C, to allow the bromine to 
convert the lead tetramethyl to lead bromide, a 
procedure similar to that used by Leighton and 


® Cf. May, Taylor, and Burton, J. Am. Chem. Soc. 63, 
249 (1941). 


vacuum 


Fic. 1. Arrangement for removing PbR,. 


Mortensen® and others for similar purposes, 
such as the determination of lead in lead tetra- 
ethyl for gasoline.'° The volatile substances 
(including carbon tetrachloride and excess of 
bromine) were then all distilled back, slowly and 
carefully to avoid bumping, at room temperature, 
into the original trap Z now cooled again to 
—196°, leaving behind just the lead bromide 
corresponding to the lead picked up in the run. 
More rapid distillation at higher termperature 
sometimes resulted in bumping of the mixture 
and loss of sample; and more prolonged treat- 
ment with bromine produced, evidently by 
reaction with acetone, an oily material very 
difficult to separate from the lead bromide. 


Analytical Method 


For the analytical determination of the lead 
bromide thus obtained, the contents of the trap 
M were transferred, after heating with 15 cc of 
very dilute HNO; (~0.3 percent), to a 125-cc 
separatory funnel, the trap being washed out 
with a further 15-cc portion of water. 

The lead in the solution was determined by the 
dithizone method using approximately the 
procedure described by Moskowitz," all reagents 
and apparatus being prepared essentially ac- 
cording to his directions. After adjustment of the 
pH of the aqueous solution to 8.5-10, using 


9 P. A. Leighton and R. A. Mortensen, J. Am. Chem. Soc. 
58, 453 (1936). 

K. Dosios and J. Pierri, Zeits. f. anal. Chemie 81, 214 
(1930). 

1S. Moskowitz, J. Ind. Hyg. Toxicol. 20, 457 (1938). 
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phenol red and concentrated NH,OH, two ex- 
tractions of the lead as lead dithizonate are 
carried out. The first extraction uses 8 cc of 
dithizone-chloroform solution and 5 cc of chloro- 
form, the water layer being further washed witha 
2-cc portion of chloroform. In the second ex- 
- traction, 4 cc of the same dithizone solution and 
5 cc of chloroform are used, followed again by 
washing. The amount of dithizone used repre- 
sents altogether approximately a_ threefold 
excess. The combined chloroform solution, 
containing the lead dithizonate together with 
excess of dithizone, is now treated with suc- 
cessive 100-cc portions of 0.5 percent aqueous 
KCN for the removal of the excess of dithizone. 
The complete removal of this reagent is indi- 
cated when a fresh portion of KCN solution 
remains colorless after shaking with the chloro- 
form layer. Complete removal was further veri- 
fied in some cases by testing the final KCN layer 
with 1 cc of chloroform and 2 drops of lead 
nitrate solution equivalent to iyg of lead. 
The lead dithizonate is next decomposed by 
treatment with about 100 cc of 1 percent HCl, 
which removes the lead and leaves the green 
dithizone in the chloroform. 

This free dithizone, equivalent to the lead 
being determined, is now titrated with stand- 
ard lead nitrate solution containing 10yug of 
Pb per cc, freshly prepared from a solution of 
1000ug/cc. The separated chloroform-dithizone 
layer is treated with 50-60 cc of fresh KCN 
solution, so that the dithizone is distributed in 
both phases, its concentration being actually 
very much higher in the aqueous KCN than in 
the chloroform layer. To this two-phase system 
is now added about three-fourths of the required 
lead, estimated after a preliminary titration of a 
similar sample or from the result of a similar run. 
The red lead dithizonate formed enters the 
chloroform and the excess dithizone remains 
practically entirely in the KCN layer. After 
removal then of the main chloroform layer, the 
KCN solution is finally titrated with small 
portions (1 cc chloroform and 1iyg of Pb) until 
the fresh chloroform portion is no longer colored 
pink after vigorous shaking with the KCN 
layer. 

The procedure was tested on known quantities 
of lead (5—50yg) in the form of lead nitrate solu- 
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tion. On 17 such tests the error ranged from 
—1.5 to —5Syg, averaging —3.2+1.5ug, with 
no clear correlation between the error and the 
amount of lead involved. With a correction of 
+3.2ug, the lead determinations should be 
expected to be correct within an uncertainty of 
about +1.5ug. As a further test, a sample of 
lead tetramethyl obtained by fractional distil- 
lation of an old sample of some material furn- 
ished by the Ethyl Corporation” was analyzed 
by the procedure described, including freezing 
out in trap L, distillation over to the bromine 
trap M, and removal of volatile substances and 
subsequent analysis of the lead bromide. For 
five aliquots of presumably 28.5yug taken, the 
result obtained, including the correction of 
3.2ug, was 26.6+1.0ug. The deficiency may have 
been caused by low purity of the lead tetra- 
methyl, but it seems to indicate that the de- 
terminations are easily within 10 percent and 
probably within 5 percent of the truth. Further 
evidence of their dependability, or at least of 
their relative consistency, is found in the con- 
sistency of the results obtained in the deter- 
mination of the half-life of the methyl radicals 
produced in the photochemical decomposition 
of acetone in the experiments described below. 

In order that the results obtained should have 
quantitative significance, the conditions of 
operation were standardized and evaluated in 
various respects. 


Duplication of Results 


Blank experiments (without illumination, or 
in the absence of a sensitive mirror) indicated 
that no lead was picked up beyond the limits of 
experimental error (~1.5ug). On the other hand, 
Table I indicates that under fixed conditions of 
illumination, and with a sufficiently long mirror, 
at a fixed distance from the irradiated zone, the 
amount of lead picked up in a given time was 
invariable within the limits of experimental 
error. Furthermore, the same table shows that 
for fixed conditions of illumination, rate of flow, 
and mirror position, the amount of lead picked 
up was proportional to the time of illumi- 
nation. 


12 We are indebted to Dr. G. Calingaert of the Ethyl 
Corporation for the gift of this material. 
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TABLE I. Effect of mirror length and number of mirrors on 
lead pick-up; photolysis at 45°C. 


Mirror Time, 
min. 


Mirror length 
mm 


2 mirrors, >32 
2 mirrors, >32 


* Mirror distance is the distance between nearest edges of irradiated 
zone and mirror. 


Mirror Length 


From Table I it is also seen that for the rather 
short, fixed time (5 min.) used in the experiments 
described below, the amount of lead picked up 
from mirrors =32 mm in length was not greater 
than that from the shorter mirrors (19.5 and 24 
mm). Similar experiments with mirrors at a 
distance of 18.5 cm from the irradiated zone 
showed, in agreement with these figures, that a 
mirror about 17-20 mm long was sufficient to 
stop all radicals. The actual mirrors used in the 
acetone experiments here reported were all 
between 32 and 38 mm long, easily sufficient to 
stop all mirror-active particles passing over such 
a surface. 


Stray Pick-Up 


An interesting aspect of this and other work 
using this method is the possibility that a 
considerable contribution to the lead pick-up 
during a run is from the stray or ‘‘shadow”’ 
mirror left in the length of the reaction tube 
during successive mirror depositions. Conse- 
quently, an experiment was run in which a 
mirror was driven up the tube toward the acetone 
inlet before the irradiated zone. Under these 
conditions the only lead that could be picked 
up during a run would have to come from the 
“shadow mirror.”” It was established that the 
amount of lead so picked up was negligible 
(i.e., <~1.5yg). 


Effect of Second Mirror 


A corollary experiment is to lay down two 
heavy standard mirrors, one near the irradiated 
zone and the other some distance away. As 
indicated above, it had already been established 
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that the first mirror was able to stop all mirror- 
active particles. Nevertheless when an additional 
standard mirror was placed further down the 
tube at a sufficient distance, increased pick-up 
was detected in the 45° photolysis, as shown in 
runs 28 and 30, in Table I. The explanation is 
simple. Evidence adduced by Spence and Wild 
and Herr and Noyes," as well as in this work, 
indicates that at least part of the reaction 
proceeds via 


h 
(CH,),CO—>CH;+CH,CO. (1) 


The acetyl radical is fairly stable below 60° 
and breaks up slowly either in the gas phase or on 
the wall. Apparently it does not itself react 
with metal mirrors."* Given time, however, it 
decomposes (i.e., as it passes down the reaction 
tube) in the over-all reaction: 


CH;CO—CH;+CO (2) 


and the additional free methyl radicals so 
produced can pick up lead in detectable quanti- 
ties from the second mirror. A similar experiment 
with photolysis at 100° (in connection with fur- 
ther work under way on the decomposition of 
acetone) indicates, as expected, no such addi- 
tional pick-up of lead at the second mirror. 


Half-Life Experiments 


It follows, of course, that the so-called half- 
life of methyl radicals as determined on radicals 
obtained from the photolysis of acetone at low 
temperatures (below 60°) is illusory and repre- 
sents an average of two different effects: a 
disappearance of methyl radicals during their 
passage down the reaction tube, presumably 
giving ethane, and a production of methyl 
radicals from reaction (2). 

For this reason, it was decided to illustrate the 
sensitivity of this analytical method of measuring 
the Paneth effect by half-life experiments on 
acetone photolyzed at 45°, at which temperature 
the CH;CO radicals have considerable persist- 
ence, and at 100°, at which temperature the 
CH;CO radicals produced by reaction (1) de- 


13 R, Spence and W. Wild, J. Chem. Soc. 352 (1937). 
4 —D.S. Herr and W. A. Noyes, Jr., J. Am. Chem. Soc. 62, 
2052 (1940). 
% Cf. T. G. Pearson and R. H. Purcell, J. Chem. Soc. 
1151 (1935) for the cases of As, Sb, and Te. 


| | 
Exp. Pb 
f Ne 
e 14 19.5 0 5 25 
f 15 24 0 5 24.5 
58 >32 0 7.5 35 
63 $32 0 3 15 
28 0, 12 5 28.8 
. 30 0, 12 5 30.5 
e 
d 
e 
of 
e 
i- 
d 
or 
of 
r- 
Is 
yn 
ve 
of 
in 
or 
ed 
of 
d, 
of 
or, 
he 
tal 
Ww, 
ed 
ni- 


FELDMAN, RICCI, AND BURTON 


a 


ug Pb (log scale) 
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Fic. 2. The logarithm of the active concentration of 
radicals (measured by the quantity of lead picked up) as a 
function of the time of transport. Upper curve, photolysis 
of acetone at 45°; lower curve, at 100°. 


compose via (2) presumably as rapidly as they 
are formed. 

All the individual points obtained in this way 
are listed in Table II and plotted in Fig. 2. 
In Fig. 2,2 represents the time in seconds taken 
by radicals to go from the irradiated zone to the 
mirror; 2 is calculated according to the formula 
of Paneth and Lautsch.!® The vertical axis 
represents the activity of the radicals in terms 
of the number of yg of lead picked up. The 
apparent half-life calculated from the curve 
is 1.1 (+0.1) X10- sec. for the radicals obtained 
from acetone photolyzed at 45° and 8 (+0.5) 
X10-* sec. for those produced at 100°C. These 
values are larger than any previously reported. 
The value for the half-life of methyls from the 
thermal decomposition of lead tetramethyl 
in He as a carrier gas, reported by Paneth and 
Lautsch,"® is 5.8 —8.8 X 10-* sec. at pressures of 2 
and 2.5 mm, but with no correlation with the 


pressure. Any comparison of such half-life 


values, except as determined in the same appara- 
tus and with the same technique, is probably 
unjustified, and it would seem that no absolute 
significance is to be attached to the values them- 
selves. In the present experiments for example, 
a small but systematic error in the determina- 
tion of the lead transported would have a definite 
effect on the half-life calculated. The difference 
in slopes here reported, however, (Fig. 2) indi- 
cates a real difference between the radicals 
produced at the different temperatures, the 
lower half-life value at 100° being expected, as 
explained above. 

The lower value for the number of ug of lead 


16H. Paneth and W. Lautsch, Ber. 64, 2708 (1931). 


picked up at ‘‘zero distance’ in the 45° experi- 
ments as compared to the 100° experiments is 
to be attributed, of course, to the fact that not 
all the acetyl radicals are immediately decom- 
posed at the lower temperature. No direct com- 
parison with the corresponding 100° value (lead 
pick-up) may be made, however, because the 
concentrations of acetone in the irradiated 
zone, and hence the amounts of decomposition 
per unit time, are different in the two sets of 
experiments. 

The slopes of the lines, however, are properly 
comparable for it was found that a change in 
temperature of the irradiated zone from 25° to 
100° was without effect on the weight of acetone 
transported per unit time. It may finally be 
remarked that if the points plotted in these 
graphs had been averages of a number of widely 
divergent determinations (as is usually the case 
in Paneth mirror experiments) the distinction 
between the slopes of the two lines would not 
have been justifiable. 


Steady State Concentration of Radicals 


The actual concentration of free methyl 
radicals at any point in the tube may be calcu- 
lated from the amount of lead picked up in a 
given time at that point and the quantity of gas 
passing that point in the same time. In illustra- 
tion, assuming a temperature of 300°K and a 
pressure of ~1.5 mm in the gas flowing by at 


TABLE IT. Half-life measurements; all runs 5 minutes. 


Dist. from Time of 
Temp. illumination, transport, Pb picked up, 
cm X105 sec. 


45 
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: 0 1 2 3 4 5 6 j 
30 
25 
20 
O ° 15 
| 104 0 0 24.0 
103 6 1 22.2 
102 8.5 1 20.8 
101 12.2 2 20.5 
| 100 15. 3 22.2 
107 17. 3 19.2 
99 18. 3 18.4 
| 98 21.5 4 19.3 
| 106 24.5 4 18.5 
| 97 25.5 4 18.2 
_ 96 29. 5 16.2 
| 105 31 5 17.2 : 
100 115 0 0 26.1 j 
| 114 6 1 24.0 
113 10 2 21.9 
112 12.5 2 21.8 
111 16.7 3 20.6 
. 110 24.2 4 17.1 
108 32 15.7 
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zero distance from the irradiated zone, the vol- 
ume of gas passing this point in 5 minutes 
is found (from the streaming velocity) to be 
2.7X10° cc: the number of moles of CH; picked 
up in that time at that point is 4X24 10-*/207 
or 4.6X10-7. The concentration of CH; is 
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therefore 1.7X10-" mole per cc, or 2.110-'% 
mole percent in the gas. 

Other examples of the usefulness of this 
method of counting free radicals will be given in 
forthcoming papers on the mechanism of the 
photolysis of propionaldehyde and of acetone. 
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Angles at the points of contact of a number of liquid-liquid-air systems were measured. A 
new method of measuring the interfacial tension is presented. The work of adhesion between 
the phases is determined. A new theory of spreading is offered. 


EASUREMENTS of contact angles at 

liquid-liquid-solid and liquid-solid-air in- 
terfaces have been made by a number of workers. 
Difficulties due to poor reproducibility of solid 
surfaces and due to the effects of adsorbed and 
condensed layers on solid phases, however, 
markedly affect the experimental results.'* 
Notwithstanding these difficulties the measure- 
ments of contact angles have been of consider- 
able value in both scientific and industrial 
investigations. 

In view of the importance of this work it 
is perhaps surprising that no one has made use 
of the corresponding phenomenon at liquid- 
liquid-air interfaces. This phenomenon can be 
easily investigated by bringing an air bubble to a 
liquid-liquid interface and measuring the con- 
tact angles of the system. The surfaces involved 
are definite and reproducible. No condensed 
layers of an extraneous phase can be present 
under the condition of the experiment. In addi- 
tion the systems involve surface energies for 
which values are known or can be determined. 
(The surface energies of solids cannot readily be 
determined.) 

Such measurements at liquid-liquid-air inter- 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

1N. K. Adam, Physics and Chemistry of p. 180. 

*F. E. Bartell and P. H. Cardwell, J. Am. Chem. Soc. 
494 (1942); 1530 (1942); 1641 (1942). 


faces can be used to determine the interfacial 
tension free from inexact assumptions concern- 
ing the value of the contact angles, and also to 
determine the values for the work of adhesion 
and the work of spreading of different phases. 
It is with this valuable and thermodynamically 
sound tool for the study of the energy relations 
at liquid-liquid interfaces that this paper will be 
concerned. 

The method used is similar to that of Taggart 
et al. applied previously only to the measure- 
ment of contact angles at liquid-solid-air inter- 
faces.* An air bubble is brought by means of a 
curved pipette to a hollow glass tube (3-mm 
inside diameter) sealed to a glass rod. The air 
bubble is then brought down to the interface 
of the two liquids and the angles of contact at 
the points where the interfaces meet are de- 
termined. This simple apparatus can be easily 
refined in some such manner as indicated by 
Bartell? to allow for convenient measurement 
on either side of the interface. Photographs were 
made of the air bubble, the contact angles being 
determined by direct measurement with a pro- 
tractor either upon the photograph itself or 
upon the ground glass focusing screen. The con- 
tact angles are sharply defined and reproducible, 
and are the same on approaching the interface 


3 Taggart, Taylor, and Ince, Trans. A.I.M.E. 87, 285 
(1930). 
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from either direction. Measurements made at 
different times agree with each other within the 
experimental precision of the measurement. 

Equilibrium at the interfaces of the systems 
studied was obtained quickly and did not change 
with time. This is not unexpected as at the im- 
mediate vicinity of the surfaces of discontinuity 
the phases are quickly saturated with each other. 
It should be noted also that as predicted by the 
Gibbs relation all the interfaces (liquid 1-air, 
liquid 2-air, and liquid 1-liquid 2) will be richer 
in the material for which —0o/dc is greatest. 

The only factors that can affect the angles are 
those which change the relations of the three 
interfacial. tensions. As distortion of the air 
bubble at the interface can have no effect on the 
value of the interfacial tensions it is not possible 
for distortion to change the values of the angles. 
Therefore the size of the bubble and the manner 
in which it is held against the interface do not 
affect the values obtained. 

Consider a system where the air bubble makes 
a definite contact angle with the liquid-liquid 
interface. On bringing an air bubble to a water- 
aniline interface equilibrium is established defi- 
nitely and suddenly. When the points of contact 
of the three interfaces are in equilibrium the 
following relations may be written.‘ (See Fig. 1.) 


Fw, an A COS (1) 
Ow, A SIN O—Gan,a SiN d=o,=0. = (2) 
From (1) and (2) we get 
dean, w= —ow, aL sin (8+¢)/cos ¢ 
cos 2¢ 
—[cos (@—¢)/cos ¢ ]dow, 4 -| a. (3) 
cos 


Aniline 


Fic. 1. 


4 ¢=interfacial tension; w= water; an=aniline; A =air; 
o¢,=0=tension at the points of contact at equilibrium. 
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Fic. 2. Water-air-aniline; temperature 22°C. 


Angle @ shall always be the angle measured 
through the top liquid, angle ¢ the angle through 
the bottom liquid, and angle 8 the angle through 
the gas phase. 

Equation (3) will be discussed further in 
another paper. 

The experimental procedure employed was 
to measure angle @ and then to calculate angles 
¢@ and 8 by the use of Eq. (2). Angles ¢ and 8 
should not be measured directly as the image of 
the aniline-air interface is distorted by the curva- 
ture of the interfaces and by the difference in the 
refractive index of the different phases. 

Reference to the photographs (Figs. 2 and 3) 
of the aniline-water-air systems’ show that 
6=131.1+3.0°.§ 

An air bubble when released below the inter- 
face in the aniline phase will rise to the interface 
and (if not too large) be held there showing the 
same relations of the angles within the experi- 
mental precision of the measurement. 

It was observed that over a range of tempera- 
ture from 20° to 50° there was no measurable 
change in the values of the angles.” 8 

Since a change in temperature over this range 
had no measurable effect on the values of the 
angles, we may therefore use the angles observed 
to determine the interfacial energies over the 
range of temperature for which the angles are 
constant. 


5 The aniline was practically colorless. C. P. Bakers 


analyzed. 
* Result of twenty measurements. Mean=131.1° a.d. 


7 Investigations are now being made with systems in the 
region of critical mixing. 

87. W. Wark and H. B. Cox [Trans. A.I.M.E. 112, 189 
(1934) ] report no measurable change in contact angles at 
liquid-solid-air interfaces for AT of 10 degrees. 
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Fic. 3. Water-air-aniline; temperature 33°C. 


By Eq. (2) using the surface tension values of 
Reynolds,’ 


0=ow,4 SIN O—Gan, a SiN ¢, 


46.4 
sin ¢=—— sin 6 
42 


o=56°, 
and by Eq. (1) 


0=Gan, A COS 4 COS O+Gan, 


Fan, = 6.90.6 erg/cm?. 


This value although precise only to about 
9 percent differs from the result predicted by 
Antonow’s rule (which rests on the assumption 
that @=180° and ¢=0°) by about 64 percent 
and from the value of Reynolds (who neglects 
consideration of the angle that the interface 
makes with the third phase) by about 44 
percent. 

In the case of the benzene-water-air (also 
ether-water-air) system ‘the air bubble when 
brought down to the interface through the 
benzene is not held at the interface. If pressed 
against the interface the bubble will tend to 
distort and force the interface down. On the 
other hand when the bubble is brought up to 
the interface from the water layer the benzene 
enfolds the air bubble, i.e., the air bubble goes 
into the benzene layer immediately. If an air 
bubble is brought below the interface into the 
water layer and a drop of benzene placed in 
contact with the air bubble, the benzene spreads 
out over the gas phase replacing the water in 


*W. C. Reynolds, Trans. Chem. Soc. 119, 461 (1921). 
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contact with the air so that the angle of the 
benzene-water interface against the benzene- 
air interface is 0° and that of the water-air 
interface against the benzene-air and _ the 
benzene-water interfaces is 180°. This behavior 
can be interpreted to mean that the air bubble 
is more easily “‘wetted’”’ by the benzene than by 
water. 

Considering the contact angles as the angles 
that the tangents to the curve (benzene-air and 
water-air) make with the liquid-liquid interface 
we have at the points of contact (Fig. 4) @=0° 
and ¢=180° and by Eq. (1) 


cos cos 180° 
(4) 


The same relationship holds for the ether- 
water-air system. These results are those pre- 
dicted by Antonow’s rule. Antonow’s rule is an 
empirical equation found to hold true for a 
number of systems. It can be seen that Antonow’s 
rule is a special case resulting when 6=0° and 
¢= 180° and Eq. (1) reduces to Eq. (4). 

For the nitrobenzene-water (Fig. 5), carbon 
tetrachloride-water, chloroform-water, and ma- 
lonic ester-water systems with air we have the 
condition that @=180° and ¢=0°. Applying 
Eq. (1) we see that in all these cases (R= organic 
liquid) 


TB, w= Tw, A— OB, A+ 


and Antonow’s rule will also hold. 

In all these systems an air bubble (if not too 
large) when released below the interface in 
the organic liquid will rise to the interface 
and be held there showing the same relations 
as above within the experimental precision of 
the measurement. 


WORK OF ADHESION OF THE GAS PHASE AT THE 
INTERFACE OF TWO LIQUIDS 


The work of adhesion between a gas and a 
liquid, both in contact with another liquid, at 
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constant temperature and pressure is numerically 
equal to the free energy required to separate the 
gas phase from the liquid phase. This value can 
now be determined. Let 1 represent the top 
liquid; 2 the bottom liquid; o; the initial 
tension at the point of contact; ¢,=0, the equi- 
librium tension at the points of contact. 

As the air bubble is brought down to the 
interface, at the moment of contact 6;=0°: 
i= 180°; 180°. 


(a) 


o;=02,4+01,2 COS $;+01,4 COs 


At equilibrium 
(b) 
Solving for 


COS 4 COS de. 


o1,2= —[o1,4 COS 0.402, 4 cos¢, | 
and substituting in (a) we get 
COS 4 COS ge. (5) 
Since dG= VdP —SdT — Ado 
dGr, p= — Ado, 


AGr, P= —A (o.— 
and since 


o.-=0, 
(AG/A)7, p=0;= —W/A = —work of adhesion. 
In the benzene-water-air system and the 
ether-water-air system 0,=0°; ¢,=180°; and 


(AG/A)r, p=0. The work of adhesion of the 
air bubble to the water in the presence of the 


Fic. 5. Water-air-nitrobenzene; temperature 35°C. 
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organic liquid is zero and the air bubble will not 
be held by the water. 

It is also necessary to consider the energy 
holding the gas phase to liquid 1. This can be 
calculated in the following manner. 

If the air bubble approaches the liquid-liquid 
interface from below at the moment of contact 
o;=0°; 180°; 180°. 


(c) =01,4+02,4 Cos Bi +01,2 COs 4;, 

At equilibrium 

(d) 

Solving for 


—01,2=02,4 COS +01, 4 COS 
and substituting in (c) 
cos de+o1,4 cos @,. (6) 


In the benzene-water-air and the ether-water- 
air systems ¢,=180°; @.=0°; 6,=180° and 
from (6) 


(AG/A)r, p will be negative and the gas phase 
will be held by the organic liquids (Table I). 

In the carbon tetrachloride-water-air, nitro- 
benzene-water-air, chloroform-water-air and ma- 
lonic ester-water-air systems 0,=180°; ¢.=0°; 
B.= 180°. Applying Eq. (5) 


and from Eq. (6) 
o;=0. 


The air bubble will thus be held by the bottom 
liquid but not by the top liquid (Table I). 
In the aniline-water-air system by Eq. (5) 


(AG/A)r, p= —11.1 ergs/cm? 
and by Eq. (6) 
(AG/A)7, p= —2.7 ergs/cm?. 


These values are equal to the free energy re- 
quired to separate the gas phase from the 
aniline phase (saturated with water) and from 
the water phase (saturated with aniline), 
respectively (Table I). 
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ot TABLE I.* 
Work of adhesion Work of adhesion Work of adhesion 
c Liquid 1 Liquid 2 _ of liquid 1. _ of liquid 2. _in air o' 
(saturated (saturated (in contact with (in contact with liquid 1 to 
be with with liquid 2) to air liquid 1) to air liquid 2 Temperature 
liquid 2) liquid 1) =(-—AG/A)7, p =(—AG/A) 7, p =(-—AG/A) 7, p 
‘id Benzene Water 68.8 0.0 57.6 19 
t Ether Water 21.2 0.0 35.0 18 
= Water Aniline 2. 11.1 81.7 26 
Water Chloroform 0. 66.8 52.8 18 
Water Carbon 0. 87.0 53.4 17 
tetrachloride 
Water Nitrobenzene 0. 49.4 86.4 18 


* Units of Tables I and II =ergs/cm?. 


Values in Tables I and II based on surface tension data of Reynolds (see reference 9). 


WORK OF ADHESION OF TWO-LIQUID PHASES 
AT A LIQUID-GAS INTERFACE 


The work of adhesion between two liquids in 
contact with a gas phase at constant temperature 
and pressure is numerically equal to the free 
energy per unit area required to pull the liquids 


—202,4 (again twice the surface energy of the 
organic phase). 

For the aniline-water-air system from Eq. (7) 
the work of adhesion of the two phases in the 
presence of air is found to be 81.7 ergs/cm? 
(Table I). 


(6) apart. This can also be calculated. In the manner indicated above the value for 
When a drop of the top liquid first touches the the work of adhesion between any two liquids 
ter- surface of the bottom liquid 8;=0°; @;=180°; in the presence of a gas phase can be de- 
and ¢;= 180°; and termined. The values thus obtained (Table I) 
do not agree with those given by Harkins"! 
(e) 9;=02,1+01, 4 COS 0;+02,1 COS $i, which are determined by measurements of the 
surface energies of the pure liquids, the inter- 
facial tension between the saturated phases, and 
_ At equilibrium the application of the Dupré equation. If the 
surface tensions of each phase saturated with 
0°: each other were determined, and the interfacial 
: Solving for tension could be independently determined 
accurately, the application of the Dupré equa- 
72,1= — 01,4 COS O,—a2, 4 COS de, tion should then give the values as above calcu- 
and substituting in (e) we get! lated (Table 1). 
CONDITIONS WHICH WILL 
DETERMINE SPREADING 
—[o1,4 Cos cos (7) The theory of spreading which is commonly 
) In the benzene-water-air and the ether-water- 2ccepted states that spreading will occur if the 
air systems 0,=0°; ¢-=180°; B.=180°; and work of adhesion between two phases is greater 
from (7) than the work of cohesion of the phase which is 
spreading.” 8 From Eq. (7) we can see that for 
(AG/A)r, p= —201,4. the organic liquid-water-air systems where the 
In the nitrobenzene-water-air, carbon tetra- Contact angles are 0° and 180° the work of 
chloride-water-air, chloroform-water-air, and ma- @dhesion will be equal to twice the surface 
y ya lonic ester-water-air systems 6,=180°; ¢-=0°; tension of the liquid which has the lower surface 
from Be= 180°. (AG/A)r, p will therefore be equal to, wy Harkins, J. Am. Chem. Soc., 356 (1917). 
line), —— 2 W. D. Harkins, 6th Colloid Symposium Monograph, 


?Compare with the Dupré equation, W=o2,4+01,4 


p. 23. 
18 W. D. Harkins, J. Chem. Phys. 9, 552 (1941). 
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TABLE II. 
(-AG/A) 7, p (-AG/A) p, p 
Liquid 1 Liquid 2 liquid 1 “Spreading liquid 2 “Spreading 
(saturated (saturated spreading coefficient” spreading coefficient” 
with wit on liquid 2 (Harkins) on liquid 1 (Harkins) Temperature 
liquid 2) liquid 1) in air t=20°C in air ¢ =20°C i 
Benzene Water 57.6 8.86 —11.2 —78.92 19 
Ether Water 35.0 44.95 13.8 — 66.35 18 
Water Aniline 70.6 —35.94 79.0 24.40 26 
Water Chloroform — 14.0 — 78.25 52.8 12.99 18 
Water Carbon —33.6 — 89.35 53.4 2.83 17 
tetrachloride 
Water Nitrobenzene 37.0 — 55.03 86.4 3.71 18 
tension. Now the work of cohesion of a liquid SUMMARY 


in air is twice the surface tension. Therefore in 
considering the liquid with the lower surface 
tension as spreading, the “‘spreading coefficient” 
defined as Wa-Wce (for the cases above mentioned 
where the contact angles are 0° and 180°) will 
not be greater than zero. When the liquid with 
the greater surface tension is considered in 
spreading the “spreading coefficient’”’ as defined 
will always be negative. 

This ‘‘peculiar restriction” in the instances 
where the empirical relation expressed by An- 
tonow actually does hold true is realized by the 
proponents of that theory, in answer to which 
they express great doubt that the Antonow 
relation is a valid one. 

Consideration of what happens when a liquid 
spreads makes it evident that to explain spread- 
ing in a specific gas phase, we need only compare 
the work of adhesion of the gas to the liquid on 
which the spreading is to occur with the work of 
adhesion of the two liquids in the presence of that 
same gas phase. If the free energy change at 
constant temperature and pressure is negative 
the action will be spontaneous, and spreading will 
occur in that direction which will give a negative 
free energy change. A direct comparison of the 
values of the free energy change for the systems 
involved will then give a thermodynamically 
exact measure of the tendency and direction in 
which spreading will occur (Table II). 

This theory is equally well applied to the 
phenomena at liquid-solid-air and liquid-liquid- 
solid interfaces. | 


1. A method is presented whereby it is 
possible to determine the interfacial tension 
between two liquids from the surface tensions. 

2. Antonow’s rule is shown to be a special 
case of the equation 


0=01, cos O+024 cos ¢. 


3. The Dupré equation is indicated as being a 
special case of the equation 


COS 4 COS J. 


4. The work of adhesion of a gas to each of the 
two liquids in contact with each other is calcu- 
lated by the equations 


+01, 4 COS 4 COS 
and 
COS +01, 4 COS 


5. The work of adhesion of two liquids in 
contact with air is determined by the equation 


A COS 4 COS |. 
6. A new theory of spreading is offered. 
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The method of the grand partition function may be 
used to calculate distribution functions F,(z, {m}) pro- 
portional to the probability that m molecules in a system 
of fugacity z, and fixed temperature, occupy the position 
of their coordinates symbolized by {n}. The method 
makes use of the distribution functions F,,(0, {n}) at zero 
fugacity. The distribution functions may be written 
F,(z, {n}) =exp [— W,.(z, {n})/kT], in which W,(z, {n}) 
is the potential of average force of m molecules at the 
fugacity z, which becomes equal to the ordinary potential 
energy at zero fugacity. The equations may be generalized 
to permit the calculation of the distribution functions at 


any fugacity assuming a knowledge of them at any other 
fugacity. Using methods previously employed for imper- 
fect gases, the pressure, and also the density in molecules 
per unit volume, may be developed in a power series of 
difference of fugacity around any arbitrary fugacity. The 
coefficients of these developments are calculable at all 
fugacities by the same equations always employing the 
potentials of average force at the fugacity around which the 
development is made. The power series obtained represent 
functions which are regular on the real positive axis of 
fugacity except at the points characteristic of the phase 

transitions in the system. 


INTRODUCTION 


E consider only systems of infinite extent 

in the volume space, the temperatures of 

which are fixed at 7. The systems are one 

component systems, containing only one kind of 

molecule. Instead of using pressure P or density 

p for the independent variable of the system we 

shall find it convenient to use the fugacity, 

normalized in density units, for which the symbol 

zor y will be used. The definition of this fugacity 
may be written: 


s=lim [po exp (u—yo)/RT], (1) 


in which y» is the chemical potential of the 
system, and yo is the chemical potential at the 
density po. The density is expressed in molecules 
per unit volume, p=lim (N/V). According to 
this definition of z, as the-fugacity is reduced to 
such a value that the system becomes a perfect 
gas the fugacity and density become equal: 


lim | =1. (2) 


z=0 p 
The individual molecule requires f coordinates 
for the description of its configuration, and three 
of these may always be chosen as the three 
cartesian coordinates, x, y, z, of the center of 
mass. The f—3 internal coordinates are either 
bounded between limits such as 0 and 2z, or the 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 


internal energy of the molecule is assumed to 
approach infinity as the value of the coordinate 
becomes large. The coordinates of the molecule 7 
will be symbolized by 


(4) Vix Qriy is (3) 
and the volume element, calculated in the f 
dimensional cartesian space by 
x 
q 


in which J is the Jacobian. The coordinates of 
n molecules will be indicated by 


{n} =(1), (2), (i), (n), (5) 
and the volume element by 
d{n} =d(1)d(2)---d(i)---d(n). (6) 


The distribution function F,(z, {n}), a func- 
tion of the coordinates of ” molecules, isdefined 
by the statement that in an infinite system of 
fugacity z the probability that molecules be 
at the position symbolized by {n}, and in the 
element of configuration space d{n} is propor- 
tional to F,(z, {n})d{n}: 


F,(z, {n}) is proportional to probability 
of the configuration {n}, (7) 


and the normalization is so chosen that 


im | [fo fre in})ain} |=1; (8) 
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the average value of the function F, integrated 
over the internal coordinates of the molecules is 
unity in the 3n-dimensional cartesian space of 
the coordinates of the centers of mass. 

Functions W,(z, {n}) of the coordinates of n 
molecules in a system of fugacity z may be 
defined by the equation 


—W,(z, {n})/kT=In F,(2, {n}). (9) 


These functions are actually the potentials of 
the average forces between m molecules in an 
infinite system of the fugacity z. 

At zero fugacity the density is zero, and the 
n molecules are isolated in space. In a classical 


_ system the distribution function at zero density 


is proportional to the exponent of minus the 
potential energy divided by kT. We see, then, 
that for a classical system 


W,(0, {n})=U,{n} (classical). (10) 


If the deviations due to quantization from the 
classical equations are large the function 
W,(0, {n}) is generally a function of the temper- 
ature, but will still have such properties of the 
potential energy as will be made use of in this 
article. For convenience we shall refer to W, at 
zero fugacity as the potential energy, but the 
use of the equations of this article is in no way 
restricted to classical systems.'! For a quantum 
mechanical system the function W,(0, {n}) must 
be calculated by means of a Slatersum. An 
alternative and probably more satisfactory 
method of handling the case of a system of 
molecules having internal quantized degrees of 
freedom would be to replace the internal co- 
ordinates by the corresponding set of quantum 
numbers and replacing integration in the equa- 
tions of this article by summation over the 
discrete values of the internal quantum numbers. 
In general the functions of the coordinates of 
the centers of mass may be handled as classical, 
but the functions would be different for different 
values of the internal quantum numbers. 

In the succeeding section we shall construct 
an equation relating the distribution functions 
at two different fugacities, z and y+z. 


1 See, for instance, the analagous treatment of quantum 
mechanical systems by B. Kahn and G. E. Uhlenbeck, 
Physica 5, 237 (1938), and by J. de Boer, Dissertation, 
Amsterdam (1940). 


GENERALIZATION OF THE GRAND PARTITION 
FUNCTION 


For reasons of brevity we introduce the symbol 
o(z) =P(z)/kT, (11) 


with P(z) the pressure. 

The probability P,(V,z, {n})d{n} that in an 
infinite system of fugacity z, there will be within 
the finite volume V exactly m and no more 
molecules, located at the position indicated by 
the symbol {n} (for which all » centers of mass 
are within V), and within the element d{n} of 
the configuration space, is known, from the 
theory of the grand partition function, to be 
given by 


P,(V,2, {n})=[exp — V¢(z) J2"F,(0, {n}). (12) 


The probability that there be exactly 
molecules within V, without reference to their 
positions, is obtained by dividing the integral 
of this function over the volume V by 1!, since 
integration covers n! equivalent positions differ- 
ing only in permutations of the identical mole- 
cules. The probability of having no molecules in 
V is exp — V¢(z). The sum of these probabilities 
over all values, 0 to ©, is unity. It follows that 


exp f f f 


which is one of the usual forms of the grand 
partition function. 

The volume V must be macroscopically large 
so that the influence of surface effects is negli- 
gible. It is to be noted that the normalization of 
Eq. (8) applies only to the limit of infinite 
volume, and the integral of F, over the finite 
volume V is not, in general, V”. 

For finite volumes V and real systems with 
short range repulsive forces the sum in (13) 
always converges, since the repulsive energy at 
short distances assures us that for sufficiently 
large values of m the integrand F, approaches 
zero over the whole range of integration. Since 
the series is regular to infinite values of the 
fugacity, singularities in the pressure, indicative 
of phase transitions, appear only by going to the 
limit of infinite volume.! 
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The probability that in an infinite system of 
fugacity z, there be at the m points {n}, in the 
element d{n}, one molecule at each point, 
irrespective of the positions of other molecules 
in the system, is [p(z) ]"F,(z, {n})d{n}. This 
probability will be equal to the sum of proba- 
bilities constructed in the following way. We 
enclose the 1 points by a singly connected volume 
V. The probability that there be a molecule at 
each of the points, {n}, and exactly m other 
molecules somewhere in V will be the integral 
of Paim(V,2, {{m}+{m}})d{n+m} over all 
positions {m} within V divided by m!. The 
probability desired is the sum of these proba- 
bilities over all values of m: 


Co(2) "Fa(z, 


x Paim(V, 2, {n-+m})d{m}. (14) 


In order to obtain somewhat simpler appearing 
equations later we introduce 


) n 
F,(z, {n}), (15) 


and note that in view of Eq. (2), 


Q,(0, {n}) = F,(0, {n}). (15’) 


Introducing (15’) in (14), and using (15) and 
(12), one finds, with the completely arbitrary 
introduction of exp V ¢(0) =1, since ¢(0) =0, 


[exp {n}) —Texp ¥¢(0)] 


m=0 


x dim}. (16) 


The convention 


Qo(z)=1 (16’) 


is used. With (16’) Eq. (13) is seen to be a 
special case of (16). 

For all finite values of V, and for real systems 
of molecules with short range forces and repulsion 


at short distances, the sum of (16) converges for 
the same reason discussed in connection with 
(13). Furthermore, for finite values of V, the 
functions Q,(z, {n}) for sufficiently large values 
of m are seen to approach zero in value for all 
values of the argument {n} within the volume V. 
We form the permanently convergent sum 


m 


[exp V¢(z) ] 


I,{n} = 


m=0 m! 


Into this the expression (16) for Qnaimn(z, {n-++m}) 
is inserted: 


e wo gmtk(_)m 


I,{a} ———[exp Vo(0)] 


m=0k=0 mik 


It is seen that for any given value of m+k the 
sum of the coefficients is (1—1)"** so that only 
the term m=0, k=0, remains. It follows that 
T,,{n} =[exp V¢(0) ]Q,(0, {n}). One may write 


[exp V¢(0) ]Q,(0, {n} 


—Lexp Vo(2)] 
m. 


x f f f {n-++m})d{m}. (18) 


Let us now form [exp V¢(y+z) ]O,(y+z, {n}), 
for the fugacity y+<2, by inserting y+<z instead of 
z in the right-hand side of Eq. (16). Into this ex- 
pression Eq. (18) for [exp V¢(0) JQn4m(0, {n-+m}) 
may be inserted. The resulting expression gives 
an explicit equation for the functions at the 
fugacity y+z in terms of a double summation 
over numbers m and k of integrals of the func- 
tions at the fugacity z of n+m+k molecules. 
For constant value of m+k the coefficient is 
the sum of (y+z)"(—z)*/m!k! over all values m 
with m+k constant. The value of this sum is 
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y™tk/(m+k)!. One obtains 
[exp {n}) 


= —Lexp Vols)] 


m=0M. 


Equation (19) might be regarded as a general- 
ization of the equation of the grand partition 
function. The special case z= 0 leads to Eq. (16). 
The case y=—z leads to Eq. (18). The case 
n=0, with the convention (17) leads to an 
equation for the pressure P=kTy, Eq. (11). In 
the usual formulation of the theory of the grand 
partition function the potential energy U,{n} of 
n isolated molecules plays a singular role. In 
the formulation of Eq. (19) the fugacity zero, 
for which W,=U,, Eq. (10), is not singled out 
as exceptional, but the potentials of average 
force at any two different fugacities are 
connected. 

For finite values of V, and for real systems 
with short range forces and repulsion on close 
approach, the sum of Eq. (19) is permanently 
convergent for all values of y. There are, conse- 
quently, no singularities in the functions g or Q 
if these are calculated by using finite V values. 
There will, however, be enormously sharp breaks 
in the plot of these functions against y if macro- 
scopic volumes are used. The equation implicitly 
assumes that the limit of infinite volume is used. 

Equation (19) is entirely unsuitable, without 
transformation, for numerical solution, and for 
the same reason that the sharp breaks in the 
functions of y occur at large values of V, namely 
because for finite y and macroscopic values of V 
numerical convergence only sets in at enormously 
high terms: the largest term in the sum is that 
for which m is about 10”. In the next section 
we shall show one method, of limited value, by 
which numerical computations may be attempted. 


THE DEVELOPMENT OF Q, AS A SUM 
OF SEPARABLE FUNCTIONS 


The equations of. the last section have been 
derived almost without reference to the form of 
the functions Q,(z, {n}) for different values of n. 
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It was only necessary to assume the absolute 
convergence of the various sums appearing, 
which was assured by the assumption of a 
repulsive energy of sufficient magnitude. 

In view of their definition in terms of proba- 
bilities, all the functions Q,(z, {n}) must be 
positive. It is seen from (19) that if the Q,(0, {n}) 
are all positive, then the Q,(z, {n}) will all be 
positive for real positive values of the fugacity. 
The functions at zero fugacity are generally 
regarded as defining the system, being themselves 
defined in terms of the potential energy U,{n} 
by (10), (9), and (15’). As long as the potential 
energy is real the functions at zero fugacity are 
positive. 

Most methods of attempting to obtain the 
properties of systems in the liquid state make 
implicit use of the assumption that the potentials 
of average force, W,(z{n}), and the distribution 
functions Q,, at the fugacity z can undergo 
arbitrary small variations away from the equi- 
librium values, without leading to geometrical 
inconsistencies in their definitions as probabilities. 

Some light may be shed on whether this is 
justifiable by a consideration of Eq. (19). 
Arbitrary functions Q,(z, {n}) could only repre- 
sent a conceivable geometrical configuration of 
molecules at the fugacity z if the functions 
themselves were positive, and if the functions 
Q,(0, {n}), calculated from (19) with y=—z 
were also positive. This condition is necessary 
and sufficient. 

Now it can be independently shown that an 
arbitrary set of functions F,,{n} obeying condi- 
tions which will be set forth in the latter part of 
this section (namely, Eqs. (25), (26), and (27)) 
represents geometrically possible distribution 
functions of any finite numbers of molecules in 
an infinite system, provided the total density in 
the system is less than some value po, but are 
geometrically inconceivable as distribution func- 
tions if the density is assumed to exceed this 
value. The proof follows somewhat later. 

_ The experimental functions Q,(z, {n}) in a real 
system obey Eq. (19) leading to functions 
Q,(0, {n}) of the real system, which are, of 
course, positive. Now suppose that we were to 
use Eq. (19) to calculate certain hypothetical 
functions, Q,(—e, {n}), e>0, for which p would 
have the negative value, —a. As long as € 


; 
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remains small (but not necessarily infinitesimal), 
the functions Q,(—e, {n}) will remain positive; 
indeed, this will presumably be the case at least 
up to values of ¢ equal to the fugacity at which 
condensation occurs if the temperature is below 
the critical point. The functions Q,(z, {n}) of 
this system would then represent the distribution 
functions in another system for which the 
fugacity is z+e, and for which the distribution 
functions at zero fugacity were Q(—e), and 
which would have a higher density than the 
first. It follows that the distribution functions 
of a real system are permissible geometrical 
representations of the distribution of molecules 
at a higher density than the real density of the 
system. 

If we now only assume that the limiting 
density po, characteristic of the distribution 
functions, is a continuous function of the form 
of the functions F,{n}, then an arbitrary 
infinitesimal variation in the form of the actual 
distribution functions of a real system will not 
lead to functions which are geometrically 
unrealizable as probability functions for the 
distribution of molecules at the density of the 
system. 

It remains necessary to prove that there is one 
limiting density po for a given assumed set of 
distribution functions, F,{n}. 

Such an assumed set of distribution functions 
must first obey the conditions imposed on any 
probability functions, that 


in 


=F,1{{n}—()}, (20) 


in which the division of the integral by V is 
required by the normalization of Eq. (8). The 
set of functions will obey (20) if 


jim, Pili), (21) 


and the approach F,4; to the product F,F, is 
sufficiently rapid as the distance r,; between the 
molecule 7 and the nearest member of the set 
{n} increases. This condition will be satisfied 
for finite values of if the functions are con- 
structed according to (25) and (26), as discussed 
later. If (27) is also satisfied we have a self- 
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consistent scheme of constructing functions 
F,{n}, up to any finite value of m, which repre- 
sent possible geometrical distributions of mn 
molecules in an infinite space, and therefore at 
zero density. 

We now ask whether such a set of functions, 
defined up to some given n, can represent the 
distribution of a set of m molecules in a system 
of volume V (n/V negligibly small) in which 
there are a total of N molecules with N/V=p. 
The question is equivalent to asking whether 
there exists a positive function Fy{N} of the 
coordinates {N} of all the molecules in the 
system for which 


= ff f =F, 


But if such a positive function does exist for N 
molecules there also exists one for N—k mole- 
cules, namely, that obtained by integrating 
Fy{N} over the coordinates of k molecules, and 
dividing by V*. Conversely if no such function 
exists for N, then no such function exists for a 
number of molecules greater than JN. It is also 
trivial to see that for the type of distribution 
functions occurring in nature, namely, those for 
which repulsion between molecules exists, and 
F,,{n} approaches zero if two molecules are too 
close together, the upper limit No, for which 
such a positive distribution function can be 
found, is not, in general, infinity. It follows that 
for any arbitrarily defined set of distribution 
functions obeying (20), there exists a density 
No/V=po such that the distribution functions 
represent a conceivable geometrical distribution 
of molecules in any system of density equal to 
or less than po, but not at higher densities. 

We propose to seek solutions for the set of 
Eqs. (19) for which the distribution functions 
Q,(z, {mn} obey certain restricting conditions. 
The conditions themselves will be discussed 
first, and secondly whether the limitations are 
to be expected to be applicable to real systems. 

We wish to seek only solutions for which the 
distribution functions of two distant subsets of 
molecules {n} and {m} are equal to the product 
of the distribution functions of the two subsets 
separately. If rum is used to signify the distance, 
in the three-dimensional cartesian space of the 
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centers of mass of the molecules, of the two 
nearest molecules belonging to different sets {m} 
and {n}, then we wish to assume that always, 


jim [Onsm(, {a-+m}) ] 
=Qn(2, {n})Qm(z, {m}). (22) 


We shall use the symbol {v}, to indicate the 
coordinates of a set of v molecules, all of which 
are members of the set {n}, and })-{v}, to 
indicate summation over all possible subsets of 
the set {n} including the improper subset 
{vy},={n}. Using the functions W,(z, {n}) 
defined by (9) we may formally define functions 
w,(z, {v}) by the equation 


w(z, {n}.(—)"Walz, {n},, (23) 
for which the inverse is | 
W,(z, {n})=Lo nm(z, {v} (24) 


in which the functions w, will be referred to as 
the component potentials of average force. 

Taking the logarithm of (22), using (15) and 
(9), one sees that the assumption of Eq. (22) 
is equivalent to the assumption that 


Jim v+u))]=0, (25) 


namely, that the component potentials of average 
force approach zero in value if any one of the 
molecules of the set is far distant from all the 
other molecules of the set, or more generally, 
if the set breaks naturally into two or more 
distant subsets of molecules. 

The use of Eq. (24) to represent the total 
potential of 2 molecules as a sum of terms each 
depending on the coordinates of a subset of the 
molecules alone, is usual at zero fugacity, for 
which the potential of average force is the true 
potential of the isolated molecules. For systems 
with van der Waals forces it is usually written 
with the drastic assumption that the total 
potential is the sum of the potentials of the 
single molecules w;(0, (7)), plus the sum over 
all pairs of the mutual pair potentials 
we(0, {(7)+(j)}), which in our terminology would 
be assuming that w,=0, v=3. The assumption 
of Eq. (25) is essentially only that all forces fall 
to zero between distant molecules. 
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We wish, however, to make a more restrictive 
assumption than (25), namely, that the approach 
to zero is sufficiently rapid so that the integral 


im | ff ftom —w,(2, {v})/kT—1] 


Cexp —wilz, (26) 


t=1 


remains finite even if the volume V over which 
the range of integration is extended approaches 
infinite limits. The assumption of finite values 
for (26) will not be valid for w2(0, {(z)+(j)}) if 
the molecules 7 and j are charged. In all other 
cases of physical significance it appears to be 
likely that the assumption is justified. 

We wish actually to assume that the terms w, 
of Eq. (24) converge very rapidly, and that 
terms with »v=3 are sufficiently unimportant 
that useful results may be obtained with their 
neglect. This assumption is in no way formally 
required, but any numerical computation in- 
volving the higher terms must necessarily be 
tedious, at least. Formally, however we shall 
have to assume that 


lim [w,(z, {v})]=0 
(for all configurations {v}). (27) 


The assumptions of Eqs. (25), (26), and (27), 
will presumably be accepted as_ reasonable 
assumptions for systems composed of uncharged 
molecules at such temperatures and fugacities 
for which the liquid or gaseous phases are stable. 
The validity of the assumptions in the case that 
the crystalline phase is stable requires special 
consideration. 

For a single crystal of fixed direction of the 
crystal axes, and with fixed phase of the lattice 
points along the axes, it appears to be natural to 
accept all three assumptions as valid for the 
distribution function. Such a crystal will be 
referred to as a fixed crystal. The function 
W,(z, (z)) will be triply periodic in the coordinates 
of the center of mass of the molecule 7. The 
periodicity of this function alone determines the 
long range order in the crystal. The function 
we(z, {(t)+(j)}) would also be expected to be 
triply periodic, with the same periodicity of the 


tive 
ach 
al 


(26) 


lich 
shes 
lues 
}) if 
ther 
be 


5 Ww, 
hat 
cant 
heir 
ally 

in- 

be 
hall 


(27) 


27), 
able 
ged 
ties 
ble. 
chat 
cial 


the 
tice 
to 
the 
be 
tion 
ates 
The 
the 
tion 
. be 
the 


lattice, in a rigid translation of the two centers 
of mass of the two molecules i and j together. 
(It might, however, conceivably be accidently 
invariant under such a translation.) The function 
w, determines the extent of the correlation 
between the positions of the two molecules, and 
must, in a real system, approach zero for large 
distances between the two molecules. It must 
also, in a real system, approach positively 
infinite values as the distance between the two 
molecules becomes small, assuring zero proba- 
bility for the two molecules to be at the same 
position. 

There seems, however, to be no simple way to 
insert boundary conditions into the grand 
ensemble upon which Eq. (19) is based, in such 
a way as to assure us that the solutions obtained 
would correspond to those of a fixed crystal. It 
is true that Kirkwood and Monroe? in their 
theory of fusion have obtained periodic solutions 
for w,; by a method which is not essentially 
dissimilar to that which can be based on (19). 
However, the absence of satisfactory existence 
theorems for the type of non-linear integral 
equation which they were forced to solve, and 
the necessity which they faced of making several 
assumptions, may still leave possible doubt as 
to the validity of the result. 

It might seem at first sight as though there 
were no difficulty. The solution expected from 
such an equation as (19) would correspond to 
the distribution functions for an ensemble of 
randomly oriented crystals with random phase 
of the lattice points for each orientation. Such 
a distribution function will be a linear combina- 
tion of those corresponding to the fixed crystals. 
The solution of an integral equation derived by 
operating on (19) may well be expected to 
correspond to that of a fixed crystal with fixed 
but arbitrary axes and phase, so that the general 
solution would be a linear combination of these 
with all axes and phases. The difficulty arises in 
that the equations lead directly to the potentials 
of average force, and under the assumptions 
which we have made in Eqs. (25), (26), and (27). 

The linear combination of distribution functions 
corresponding to an ensemble of fixed crystals 
will not correspond to a linear combination of 


2 John G. Kirkwood and Elizabeth Monroe, J. Chem. 
Phys. 9, 514 (1941). 
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potentials of average force, and the potentials of 
average force of such a linear combination will 
most certainly not obey the condition (22). None 
of the assumptions (25), (26), or (27) will be 
valid for the potentials of average force obtained 
from the distribution functions for an ensemble 
of crystals by Eq. (9). 

We propose to show, however, that solutions 
corresponding to one single fixed crystal would 
be expected, if a system is assumed whose 
characteristics go over in the limit to those of 
the real system. It appears therefore to be 
probable, even aside from the success of Kirk- 
wood and Monroe, that distribution functions 
corresponding to one fixed crystal, and obeying 
(25), (26), and (27), would be solutions to 
equations based on (19) and solved with the 
usual limits of accuracy in the treatment of the 
grand partition function. 

Consider a real system, characterized by its 
component potentials w,(0, {v}) at zero fugacity, 
for which the solutions of Eq. (19) at T and z 
would be expected to be functions corresponding 
to an ensemble of crystals with randomly 
oriented axes. The component potential w;(0, (7)) 
of molecule i of this system is independent of 
the coordinates of the center of mass. Consider, 
now, a second system, consisting of a finite 
number N of molecules, but otherwise differing 
from the first only in an additive triply periodic 
term w’;(0, (2)) in the component potentials of 
each of the single molecules, which has the 
periodicity of one of the fixed crystals of the 
equilibrium ensemble of the first system at T 
and z. The average value of w’; may be chosen 
as zero, and its amplitude as proportional to 
N- In N. The equilibrium configuration of this 
system will be that of a single fixed crystal, the 
orientation of which will be determined by the 
periodic external field, and the relative proba- 
bility of this orientation to a random one will be 
proportional to some power of N, approaching 
infinity as the size of the system becomes infinite. 
Similarly the equilibrium configuration of a 
finite portion of volume V of this infinite system 
will be that of a single fixed crystal. But as the 
number N of molecules in the ensemble ap- 
proaches infinity, the functions Q,(0, {n}) of 
any finite number of molecules approach the 
same value as the corresponding functions of the 
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first system without the periodic external field. 
For solutions of (19) at finite values of V only 
functions Q, with n finite enter. We must expect 
functions Q,(z, {n}) corresponding to a single 
fixed crystal to be also solutions of (19) for the 
original system, without the periodic field, as 
well as those not obeying the conditions (25), 
(26), and (27), corresponding to an ensemble of 
randomly oriented crystals. 

The apparent contradiction of postulating 
two different solutions to (19), which in the 
form z=0, consists of a set of completely solved 
equations for the different Q,(y, {n}) functions 
of different 1 values, is resolved by the fact that 
the equation must always be solved by one of a 
number of various approximation methods. It is 
assumed that as V becomes large the influence 
of the surface is immaterial, and the ‘‘true”’ 
solutions are obtained without reference to the 
surface conditions. 

Condition (27) assumes that the functions 


Q,(z, {n}) for large but finite values of are 
completely determined in their functional de- 
pendence on the coordinates {n} by a finite (and 
numerically to a good approximation by a 
small) number of functions w,(z, {v}). 

Functions Q,(z, {n}) obeying condition (22), 
and therefore also (25), may be conveniently 
written in various ways as sums of products of 
functions which are separable in the coordinates 
of subsets {v}, of the set of molecules {n}. 

Condition (26) is required in order that the 
integral over the functions of the subsets 
converge. 

For instance, we may write Q, as the sum, 
dX {k{vi}n}u, over all possible sets, {k{vi},}u, for 
various k values, of k unconnected subsets {v,}, 
of the set {n}, of products of functions of these 
subsets {v;},. By unconnected subsets we mean 
that no member occurs in more than one subset, 
but every member of the set {n} occurs in one 
of the k subsets of every set of subsets. We write 


On(s, {m}) =E Lole)/s Tels, (28) 
for which the inverse is ; 


We define the cluster integrals as 


f 


als, (30) 


which are, in the limit of infinite V, finite and independent of V for all finite values of », in view 


of (26). 
Using (30) in (28) one finds 


v21 Ny: 


1 


since the number of ways that one can arrange m numbered objects in piles, 1, piles of v objects 
each, is n!/IIn,!(v!)" with }- »n,=n. The normalization (8) of F,, with (15), shows that 


J out {m})d{n} =[ds(2) (32) 


by(z) 1, (32’) 


or 


which fixes the arbitrary constant in w;(z, (z)) so that the normalization (8) is valid. 
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We must also be able to integrate the function Qnim(z, {m-++m}) over the limited number of 
coordinates of the set {m}. It is somewhat more convenient to write our definitions in terms of 
the functions F,,(z, {n}), related to the Q,’s by the simple proportionality of Eq. (15). We define 
functions G,,(z, {n}, {wi}m) by 


Frim(Z, {n}) TT Gru(Z, {0}, twi}m)s (33) 


for which the inverse is 
{mn}, TT {{n}+{nm,},})/Fn(z, {n})]. (34) 


If the molecules of the set {yu} and those of the set {n} are all distant in the three-dimensional 
space of the centers of mass, the functions G,,(z, {n}, {u}) approach the functions g,(z, {u}) in value. 
We make use of this latter fact and write 


(my, fv} bv} n, (35) 
for which the inverse is 
=(—)’gu(2, tn}(—) "Gauls, (36) 


The functions g,, approach zero in value if any one of the molecules of the set {u} is sufficiently 
distant from all members of the set {v}. 
We may substitute (34) for G,, in (36) and obtain 


ink 
II Fn+m(2, IT F,,.(2, (37) 


or with (35) in (33) one finds 
Frym(2, {n-+m}) = F,(z, {m}) {el IT {gulZ, ADL nBrulZ fv} ny (38) 


In terms of the functions Q, using (15), this becomes 


Onsm(2, {n-+m})=Qn(z, {n}) {ws)m) 


It is seen that with Qo=1, Eq. (28) is the special case of (38’) with n=O. 
We define 


foe tv}, tu} dtu}, (39) 


analogously to the cluster integrals of Eq. (30). Using (39) in (38’) one finds 


my 20 


f f Onim(z, {n-+m})d{m} =Qn(z, {m})m! um, =m) Vb,(2)[(2)/2} 


+X tv} n)Lo(z)/2 (40) 


It is seen again that (31) is a special case of (40) when n=0. 
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DEVELOPMENT IN POWERS OF THE DIFFERENCE OF FUGACITY 


Equation (40), if inserted in (19), leads immediately to expressions for the pressure, P(y+2z) 
=kT¢(y+z), and the component potentials of average force w,(y+z, {v}) at the fugacity y+<, 
in terms of a power series in y, the coefficients of which are the functions },(z) and 6,,(z, {v}) at 
the fugacity z. 

Multiplication of Qnim by y™ in (40) may be accomplished by multiplication of each term p(z)/z 
by y. Division by m! and summation over m (Eq. (19)), removes the restriction on the m,’s. One finds 
1 
Lexp Vo(y+z) ]On(y+z, {n})=Lexp Vo(z)]Q,(z, 


my21 m,! 
+D {v} Ve(z) ]Q,(z, {n}) 


wal pai 


for all values of y smaller than the radius of convergence of the series occurring in the last expression. 
Using (15) and the Eqs. (9) and (13) defining the component potentials of average force, taking 
the logarithm of (41) and collecting the terms on both sides of the equation that are the same power 
of V and functions of the coordinates of the same molecules, one finds: 


e(y+2) = (42) 


ual 


wi(ytz, (¢))/kT+1n [p(y +2)/(y+2) [p(2)/2]+ (43) 


pai 


w(y+z, {v})/kT=w,(z, {v})/RT+ 0,,(2, tv} )Lye(z)/2}. 


pai 


(44) 


For the special case z=0, p(z)/z=1, Eq. (2), with y+z=y, and g(z)=0, Eq. (42) becomes the 
usual equation for the pressure expressed as a power series of the fugacity in an imperfect gas.’ 
Under the same conditions Eqs. (43) and (44) permit the calculation of the component potentials 
of average force in an imperfect gas in terms of a power series in the fugacity. The equations are 
essentially similar to those derived previously‘ although the form corresponding to (43) and (44) 
has not been published previously. 

As with the special case that z=0, the singularities on the real axis of the series (42) correspond 
to the singularities in the function g(y+z2) in the real system.® 

However, singularities off the real axis of y, or on the real axis with opposite sign, may cause 
divergence of the series (42) at values of y for which no singularity in ¢(y+z) exists. Equations 
(43) and (44) present feasible methods of calculation to extend the evaluation of g(y+z) beyond 
such radii of convergence. One may use (43) and (44) to evaluate w,(y’+z, {v}) for some value y’ 
inside the radius of convergence. With these functions known, the coefficients b,(y’+z) for the 
development around the new point y’+z may be evaluated. This may be repeated until a singularity 
of g on the real axis is reached. 

One might, of course, attempt to calculate around the singularity on the real axis by following 
this procedure into the complex plane. 

Some interesting equations may be obtained by differentiation of (42). We differentiate with 
respect to y, v times, and let y approach zero. One finds 


[9’e(y+2)/ dy" ],-0= 8’ 9(2)/d2” 


3 Joseph E. Mayer and S. F. Harrison, J. Chem. Phys. 6, 87 (1938). 
4 Joseph E. Mayer and Elliott Montroll, J. Chem. Phys. 9, 2 (1941). 
5 See the discussion of Max Born and Klaus Fuchs, Proc. Roy. Soc. A166, 391 (1938). 
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This may be rewritten as 


8’ o(y+z2)/dy’ = (45) 


For the particular case that y=1 this reduces to a readily obtained thermodynamic equation, 
in view of (32’) that b;=1, 


(46) 


We may equate the expression on the right of (45) to that obtained by differentiation of (42) 
without having y approach zero, obtaining 


(u—v)!y! 
For the case v=1 we obtain 


ye(y+z2)/(y+2) = ub, (2)[yp(z)/2}, (48) 


which permits the calculation of the density at the fugacity y+z from the properties of the system 
at the fugacity z. 

In the treatment of the monatomic imperfect gas*‘ it was found to be advisable to make use of 
the fact that the integrals b, could be expressed as a sum of products of what were there called 
irreducible integrals 8,. This led to an expansion in a power series in the density p as an expression 
for the pressure, the coefficients of which involved these irreducible integrals. The transformation ~ 
had several advantages. In the first place, it is convenient to have density rather than fugacity as 
independent variable. Secondly, the calculation of the cluster integrals would actually be performed 
by primary evaluation of the irreducible integrals, and it seems natural to attempt to express the 
properties of the system in terms of these integrals. 

A third result was more complex in nature. It is possible that the power series in fugacity diverge 
while those in the density remain convergent, the criterion for the phase transition (divergence of 
the fugacity series on the real axis) being either that the density series become singular, or that one 
of them (>°k6ip*) becomes unity in value. The type of transition would be ordinary first order in 
the former case, and anomalous first order in the second case. It was claimed that the disappearance 
of the phase transition with rising temperature at the critical point must be preceded by a tempera- 
ture range in which the anomalous transition is present.* The conclusion has been criticized! on 
the justifiable grounds that the physical nature of the coefficients, the irreducible integrals 6,, 
were not made use of. The criticism, in this form, may be answered, at least partially. The 6's 
may be defined as integrals over the coordinates of k+1 molecules (divided by the volume) of an 
integrand which is a function of the potential energy and of T. Especially if component potentials 
w,(0, {v}) with v=3 are assumed to exist (which was not assumed in the previous paper) there 
seems to be no reason why physical systems should not exist in which a wide variety of ratios would 
occur between the numerical values of the 8;’s for different k’s. 

As the parameter T increases a singularity on the real axis of p in the series }>kB,p* may well 
disappear into the complex plane at T=T7,, and the series remain regular (represent a regular 
function) to high values of p on the real axis above this temperature. If the series approaches a 
value less than unity when the singularity is approached on the real axis for all T=T, then the 
density po(T) at which the singularity occurs will be the density of condensation. T, will be the 
upper limit of the temperature at which condensation occurs. ; 

However, if the series approaches a value less than unity at the singularity as T approaches T,, 
then (@P/dV) r= —(RT/V?)(1—SkBxp*) will not be zero at the critical point. This state of affairs 
would represent a strange type of critical phenomenon not found in nature. In order to have a 
critical point of the van der Waals type it is necessary that the value of the function represented 


°S. F. Harrison and J. E. Mayer, J. Chem. Phys. 6, 101 (1938). 
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by >-kB:p* (approached along the real axis from small p values) at the singularity become exactly 
unity as T approaches 7,, and that the value of the derivative }°k?8;,p* become zero at this point. 
These criteria can be fulfilled, if at all, only by very special types of potential energy functions, 
and it appears hardly possible that all critical phenomena are of this type. 

The actual state of affairs in real systems seems to be that the value of >>, p* at the singularity 
increases with T to unity at some temperature 7, lower than the temperature at which the singularity 
moves into the complex plane. Above 7, an anomalous first-order transition occurs. The experi- 
mental evidence is certainly not in contradiction to this view. 

Since the critical phenomena in binary liquid mixtures appear to show the same anomalous 
first-order transition region that is evidenced in the one component gas liquid critical region,’ it is 
tempting to try the same introduction of the 6;(z)’s in the equations derived here, in the prospect 
of eventually treating two or more component systems. 

In order to make the development general, it is necessary to define the new quantities introduced, 
Bm(z) and hnm(z, {n}), in terms of the },(z) and b,,(z, {v}), in which case the development is purely 
formal and any conclusions drawn about the nature of the approach to singularities would be open 
to the criticism previously raised. However, in the case that w;(z, (z)) is invariant under translation 
in the three-dimensional space of the center of mass of the molecule 7, which is the case in the liquid 
and gas, and that w,(z, {v}) for v=2 is independent of the internal coordinates, then the functions 
Bn(z) and hnm(z, {m}) are calculable as direct integrals of relatively simple functions of the potentials 
‘ of average force and temperature. In this case the definitions of these functions are given for 6,(0) 
in Eq. (18), reference 3 and for hyn(0, {u}) in Eq. (53), reference 4. The definitions given in these 
articles assume that we alone differs from zero, but the generalization for other cases is fairly obvious. 
In the general case we may use the definitions: 


( -i+ »)! b,( 


my>0 m! om, ! 


(49) 


for which the inverse is 


1 Bm Hm 
b,(z) =— = 1) Il Le J (50) 
um 20 Mm! 
and 
byy(2, {v}) =ub,(2)h,1(z, J. (51) 


With these definitions the method used in references 3 and 4 may be readily generalized to 
apply to the Eqs. (42), (43), and (44). Equation (49) shows that m8,,(z) is the negative of the co- 


efficient of y” in 


so that 


1 1 


m21 | 


1 
-— g E (52) 


where the path of integration must enclose the origin of y, and, in order to permit the second ex- 
pression, must be just outside the curve along which the absolute value of the sum, >-ub,(z)y", 
exceeds x in value. The integrand has two poles enclosed by this path: one at y=0 with residue 


7 Louis D. Roberts and J. E. Mayer, J. Chem. Phys. 9, 852 (1941). 
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unity, and the other when the sum equals x, that is at yp determined by 
ub, (2) yo" =x, (53) 


| 


the residue of which is found by inserting the derivative of the denominator for the denominator 
to be —x[>0wb,(z)yo" }-. Comparing (53) and (48) we may set x= yp(y+z)/(y+z) and find that 


m21 


m21 


Equation (54) is valid up to values of y for which the function g(y+z) becomes singular, and 


this y is determined by 
]"=1, (55a) 


singular. 


m21 


(55b) 


For the not too trivial proof that the case (55b) as well as (55a) determines the singularity in g(y+z) 
the reader is referred elsewhere.*® 
To proceed we note that from (53) 

d In yo/dx=> pb, (2) yo" =x-!— mB, (z)x"—!. 


| 


Integrating and substituting we find 
In [ye(y+2)/(y+2) ]=In Lyo(2)/2]+ (y+z)]”. (56) 
Similarly, using (53) again, we find 
(d/dx) b,(2)yo"=x(d In yo/dx)=1— mBn(2)x™ 


| 


which can be integrated and substituted to yield 


o(y+z) = + 


pal 


m2=1 


We shall now show that Eqs. (43) and (44) may be transformed into a power series in 
ye(y+z)/(y+z), the coefficients of which are the functions hym(z, {v}). To do this the symbol x is 
again used for ye(y+z)/(y+z) and yo for yp(z)/z. From Eq. (57) it is seen that 


(9/OBm—1)2 = (9/ABm—1) = —(m—1)x"/m. 
One may, instead, use (53) and (56) to find 
(0/OBm—1)2 (9b, /PBm—1) + (9 In yo/OBm—1)2 
=D yo" — x", 
and comparing this with the previous equation, one finds 


Lye(2)/2 = Lye(y+2)/(y+2) J"/m. (58) 


* Mayer and Mayer, Statistical Mechanics (J. Wiley and Sons, 1940), p. 465. 
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If (48) and (58) are used with (51) in Eqs. (43) and (44) one obtains immediately 
wi(y+ez, (¢))/kT+1n [p(y+z)/(y+2) ]=wilz, (¢))/kT+I1n [p(2)/z] 
+ him(z, (59) 


w(y+2, {v})/kT=w,(2, {v})/RT+ hom(z, (60) 


For the special case that z=0, y+z=y, ye(y+z)/(y +z) =p(y), and p(z)/z=1, Eqs. (54) to (57) 
reduce to those derived in reference 3, and (59) and (60) reduce to equations derived in reference 4. 


SUMMARY AND CONCLUSION 


Let us, for brevity, adopt the symbol 


in which z is the fugacity, Eq. (1), (zs) =P(z)/kT is proportional to the pressure P, and the Q,’s, 
Eqs. (15) and (7), are proportional to the distribution functions at the fugacity z. At zero fugacity 
¢(0) =0 and in view of (2), and (9) and (10), 


Jx(0, V)=Ja(V) = f f f exp (—Us{n}/kT)din}, (61’) 


in which U,{n} is the potential energy of ” molecules as a function of their coordinates {n} of 
position. 
The normal use of the grand partition function leads to Eq. (13) that 


exp Vo(z)=¥ Ja(V)z"/n!, (62) 


a power series in z with positive coefficients that for real systems with molecules having short range 
repulsion is absolutely convergent to z= © for all finite values of the volume V. The first term in 
the series, Jo, is unity. 

Differentiation of (62) leads to an equation for the density p in molecules per unit volume, 


p(2) => Jn(V)2"V/n!, (63) 
n=0 n=0 


n=0 m=0 


- (64) 
> J.(V)2"/n!]? 


n=0 

which shows, in view of the positivity of the coefficients, that neither p nor dp/dz can take negative 
values on the real positive axis of z, which alone is of physical significance. 

The first result of this paper is that if the function exp [ Vg(y+z) ] is developed around the point 
Z as a power series in y the coefficients are the J,(z, V)/n! of Eq. (61) (compare Eq. (19) with =0) 
in which the distribution functions of the integrand are proportional to exp [— W,(z, {n})/kT], 
with W,(z{n}) the potential of average force at the fugacity z; and the normalization of the distribu- 
tion functions are so chosen that their average value in an infinite volume are [p(z)/z]" exp [V¢(z) ]. 


and 
|| 


59) 
60) 


57) 


61) 


city 


61’) 


| of 


(62) 


nge 
n in 


(63) 
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CONTRIBUTION TO STATISTICAL MECHANICS 


The function g(z) may be expressed as a power series in 2, the »’th coefficient of which, },(0), 
is an integral of an integrand which depends only on the potential energy of v molecules, and the 
temperature (Eqs. (29) and (30) with z=0). The quantities b,(0) depend on the volume V over 
which the integration is extended, but for finite values of v they approach finite limits as V->~, 
provided the mutual potentials of the molecules converge rapidly enough to zero at infinite distance. 
The quantities b,(0, V) could be defined using any macroscopic but finite V, and the series 

gv(z) b,(0, Vie", (65) 


used to determine the ¢g(z) defined by (62) at this value of V. The coefficients b,(0, V) of this series 
would have the value 6,(0) of the limit V= except for very large values of v. Since (62) is absolutely 
convergent for all 2’s, if V is finite, it follows that (65) determines a function which is regular except 
at the zeros of (62), and that (65) is convergent out to the absolute value of z of the zero of (62) 
nearest to the origin. This zero of course cannot occur on the positive real axis of z. 

The function gy(z) determined by (65) is then regular along the positive real axis of z, as long as 
V remains finite. If for each v we extend the range of integration of Eq. (30) to infinity, the coefficient 
approaches the limiting value b,(0), and the resulting series 

6,(0)2” (66) 


v21 


expresses the pressure P(z)=kT¢(z) of an infinite thermodynamic system, having singularities on 
the real positive axis at the values of the fugacity characteristic of the phase transitions. At ail 
these singularities for real positive finite z's, the function g(z) represented by the series (66) ap- 
proaches the same (positive finite) value as the singularity is approached from either above or below 
along the real axis. 

The slope 0¢/dz=p(z)/z may approach a higher (positive finite) value if the singularity is ap- 
proached along the positive real axis from above than if approached from below. In this case the 
transition is first order. If dp/dz approaches finite values as the singularity is approached either 
way the transition is normal, if infinite values are approached the transition is anomalous first order. 

If the slope of g, namely, p/z, approaches the same value as the singularity is approached from 
either above or below, then the singularity corresponds to a phase transition of second or higher 
order, the order depending on which derivative of ¢ undergoes a discontinuity. 

It has been shown that if the function g(y+z) is developed around the point z as a power series 
in y, and the »’th coefficient is written [p(z)/z]’b,(z), then b,(z) is calculable by the same equation 
as b,(0), with only the replacement of the potential energy appearing in the expression for the 
latter by the potential of average force at the fugacity z. 

The component potentials of average force, w,(z, {v}), Eq. (23), may also be developed as a power 
series in the fugacity around any arbitrary fugacity. The coefficients of the developments around 
different fugacities differ only that in each case the potentials of average force at the fugacity around 
which the development is made enter into the equations. In the special case of zero fugacity the 
potentials of average force are the potential energies of the molecules. 
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An elementary systematic method for obtaining the relations between thermodynamic 
coefficients is presented. This method is applicable for any desired set of independent variables, 
and is also easily adapted to obtain the derivatives of any functions of the thermodynamic 


variables. 


DIFFICULTY invariably encountered by 
students of physical chemistry and thermo- 
dynamics is to find a systematic method for 
getting the relations between the important 
thermodynamic derivatives. One can of course 
use the tables of Bridgman,' but in cases where 
other variables, such as tension and length, enter 
the thermodynamic formulae, these tables cannot 
be used without modification. 
For the usual systems of chemical thermo- 
dynamics, the combined first and second law is :* 


dE=TdS—PdvV. 


(1) 


The free energies of Helmholtz and of Lewis 


and Randall are defined as follows: 
A=E-TS, 
F=E+PV-TS, 
H=E+PV, 


(2) 
(3) 
(4) 


(5) 
(6) 
(7) 


Since dA and dF are total differentials, we 
know that they satisfy the Euler relation. From 
(5) and (6) we can therefore derive the Maxwell 
relations (8) and (9). 


(0S/0P)r = —(dV/dT)p. 


dA =—SdT—PdV, 
dF=—SdT+ VdP, 
dH =TdS+ VdP. 


(8) 
(9) 


Since TdS is the heat absorbed in a reversible 
process, TdS. Remembering the definitions 


1See Lewis and Randall, Thermodynamics, Chapter 14. 
*We shall use the Lewis and Randall nomenclature 
throughout. 


of C, and C,; namely, 
Cy=(dQ/dT)p, 
C,=(dQ/dT)y, 


(10) 
(11) 
(12) 
(13) 


it is seen that 
Cy,=T(dS/dT)p, 


C,=T(dS/dT)y. 


The relationships discussed above are of 
course the essential core of the usual mathe- 
matical discussion of thermodynamics. In order 
to find the relationships between the various 
thermodynamic derivatives, these relations must 
be kept in mind. However, we shall discuss below 
a systematic method for obtaining the other 
thermodynamic relations, and we shall see that 
the above relations are the only ones with which 
we need burden our memory. The method is 
best explained by example. 

Suppose that we wished to express (0H/0F)s 


in terms of more familiar thermodynamic 


derivatives. First set up the equation 


dH=XdF+ Yds, (14) 


where X and Y are to be determined. X is of 
course the quantity we desire. 

We next have to specify the variables that 
we shall consider as independent. A frequent 
choice is P and T. The fundamental idea of this 
scheme is to express all the differentials in 12 
in terms of dP and dT. When this is done we 
can equate the coefficients of dV and dT on the 
right- and left-hand sides and so obtain two 
simultaneous equations from which we can solve 
for X and Y. We first substitute in (14), Eqs. 
(6) and (7). 


TdS+VdP=X(—SdT+VdP)YdS, (15) 
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but 
dS =(0S/dT) pdT+(0S/dP) 


which from (9) and (12) becomes 
dS=(C,/T)dT—(dV/dT) pdP. 
Substituting (17 in (15) 


C,dT —T(dV/dT) pdP+ VdP 
= X(—SdT+ VdP) 
(18) 
Equating the coefficients of dP and dT on both 
sides we have 


(16) 


(17) 


(19) 


(20) 
Simultaneous solution of (19) and (20) gives 
X= -, 
CpV—ST(dV/dT)p 


The extension of the method that we have 
used to other thermodynamic variables is very 
straightforward. 

Suppose that we were interested in the thermo- 
dynamics of rubber where the variables 7 and L 
representing force of tension (or compression) 
and length, respectively, were important rather 
than P and V. 

The combined first and second law can be 
written 

dE=TdS+rdL. (21) 


We define the free energies of Helmholtz and 
Lewis and Randall as follows: 


AuwE-TS. 


dA = (22) 
dF = —SdT—Lar. (23) 


From (22) and (23) we get relations correspond- 
ing to the Maxwell relations. 


(0S/dL) r= —(07/8T) 1, 
(8S/dr) p= 


(24) 
(25) 
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An interesting coefficient in the case of rubber 
is (@E/dL)r. Take L and T as independent 
variables. 

dE=XdL+YdT, 


TdS+rdL = XdL+ YdT, 


T[(S/dL) rdL + (8S/8T) dT] 
+rdL=XdL+YdT. (28) 


Introducing Eq. (24) into (28) and equating 
coefficients of dL we obtain 


(0E/dL) r= —T(d7r/dT) +7. (29) 


(From this result we can see that if 7 is a 
constant times the absolute temperature, then 
(@E/dL)r=0, which was pointed out by K. H. 
Meyer. This fact has been very important in the 
development of the statistical theories of 
rubber.) 


The method described above can be general- 
ized to more than two independent variables, 
and of course can be used to treat coefficients 
such as [0(E/T)/0(1/T) ]p or other derivatives 
involving functions of the thermodynamic vari- 
ables. For example, in the theory of fluctuations 
of density the coefficient (02S/0V")z is very 
important. From Eq. (1) (0S/8V)e=P/T. Ap- 
plying the methods used above it can readily be 
shown that 


(30) 


The second term is taken as zero for gases but 
is not at all negligible in liquids. 

Though schemes have been devised which 
successfully treat the general problem of ob- 
taining the relations between thermodynamic 
derivatives,”* the above method can be recom- 
mended for simplicity and directness. It is 
readily adaptable to whatever choice of inde- 
pendent variables may be found convenient. 

I wish to thank Professor Henry Eyring who 
suggested the need for a method of solving the 
problem suitable for the classroom. 


(26) 
(27) 


2F. Lerman, J. Chem. Phys. 5, 792 (1937). 
3 P. W. Bridgman, Phys. Rev. 3, 273 (1914). 
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In Section I arguments were presented for the necessity of splitting the log &..* term in the 
general kinetic equation into a part which represents the rate constant extrapolated to a 
standard reference state of temperature, dielectric constant and concentration, and into a part 
giving the temperature dependence. In Section II the repulsive potential between reactants at 
close distances was formulated and this potential was then incorporated in the general kinetic 
equation to give the temperature dependence of the reaction rate constant. In Section III this 
equation was applied to existing kinetic data and found satisfactory. Also the repulsive forces 
between reactants calculated from the temperature and dielectric dependence of the reaction 
rate constants were compared with the intermolecular forces between the molecules of noble 
gases calculated from the second virial coefficients as evaluated by Fowler. 


I. INTRODUCTION 


AIDLER and Eyring! point out that a term 
of the form ¢/kT should be included in the 
kinetic equation to account for the non-electro- 
static effects, which can be treated as repulsions 
beginning suddenly at the distance of closest 
approach, modified to take into account van der 
Waals attractions which may be of considerable 
importance. A similar term must account for the 
temperature dependence of the rate in the general 
kinetic equation? 


log k* =log 
(1) 


In this equation In &,,.* must be composite con- 
taining a term which represents the rate constant 
extrapolated to a standard reference state of 
temperature, dielectric constant and concentra- 
tion, and a term giving the temperature de- 
pendence. Although it is generally recognized 
that the temperature and dielectric constant 
dependences of reaction rates are quite different, 
yet Eq. (1) as stated does not indicate this 
difference. The contrast between these two 
dependences can be shown as follows: If a 
reaction between ions is studied as a function of 

1K, J. Laidler and H. Eyring, Ann. New York Acad. Sci. 
39, 303 (1940). 

2 J. N. Bronsted, Zeits. f. physik. Chemie 102, 169 (1922); 
115, 337 (1925). 
asda). Christiansen, Zeits. f. physik. Chemie 113, 35 


4G. Scatchard, J. Am. Chem. Soc. 52, 52 (1930); Chem. 
Rev. 10, 229 (1932). 


the dielectric constant of the solvent at constant 
temperature, then from the slope of the curve 


obtained by plotting log ke versus 1/D, ratrn 
can be obtained, provided that the effect of the 
solvent is that due to its dielectric constant, 
i.e., that predicted by electrostatic theory. The 
values of r4+7g obtained in these cases are 
always reasonable in magnitude and positive in 
sign. This is true because the slopes of the 


log ki-o versus 1 /D curves are of the right order 
of magnitude and of the right sign, i.e., the slopes 
are negative for ions of like charge sign and 
positive for ions of unlike charge sign. ® If Eq. 
(1) were complete, one should be able, by using 
it, to obtain r4+7rz from the slopes of plots of 
log k versus 1/T. However, irrespective of the 
signs of the charges on the reactants, the slopes 
of such plots are always negative as predicted 
by the Arrhenius equation, so that, for reactants 
of unlike charge sign r4+7r, would be negative. 
The magnitude of r4+rz calculated from a 
temperature dependence will be roughly one- 
third to one-fourth that obtained from a dielec- 
tric constant dependence for ions of like charge 
sign, but in this case the sign is positive. 

From these considerations, it follows that the 
log k.* term must be composite and must em- 
body a non-electrostatic potential which is 
predominant in determining the dependence of 

5 E.S. Amis and V. K. LaMer, J. Am. Chem. Soc. 61, 901 
(1939). 


6 W. J. Svirbely and A. Schramm, J. Am. Chem. Soc. 60, 
330 (1938). 
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the rate upon the temperature. The purpose of 
this paper is the evaluation of this non-electro- 
static potential together with a comparison of 
its magnitude with the value of electrostatic 
potentials with which we are more familiar. 


Il. THE KINETIC EQUATION 


It follows from the arguments which have 
been presented above that the mutual potential 
energy between two reacting particles can- 
not be restricted to their electrostatic potential 
energy. In order to obtain the observed tempera- 
ture dependence of the rate constant a potential 
must be added which will correspond to the 
strong repulsive forces that exist for a close 
approach of the molecules, whether these are 
charged or not. 

Following the line of thought introduced by 
Bronsted we shall assume that the reaction be- 
tween two molecules A and B takes place via 
the formation of an intermediate complex X. 
Furthermore, we assume that spherical sym- 
metry about each particle exists and that the 
mutual potential of two molecules A and B, if 
they are at the distance 7, is given by an expres- 
sion of the form 


¥(r) = gi(r) + ¢2(r). (2) 


Here ¢2(r) represents the contribution of the 
electrostatic forces in the case that A and B 
are charged (valency 24 and zz, respectively). 
Thus, under the usual assumptions, we will have 


g2(r) /(Dr(1+x«a)), (3) 


where «x is the parameter of Hiickel and Debye’s 
theory and a the distance of closest approach. 
As to the origin and form of the potential 
¢i(r) nothing need be assumed for our present 
purpose, except the following two conditions: 


1. that g; is independent of the dielectric con- 
stant D, and 

2. that g: has a positive value for those dis- 
tances r at which formation of X occurs. 


We shall call g:(r) the non-electrostatic potential 
though it may be traceable, in ultimate analysis, 
to forces of electric origin. However, in the latter 
Case, it would have to be assumed that they are 
independent of the dielectric constant. There is 


no contradiction in this assumption since, for 
the short distances involved, the notion of the 
dielectric constant (of the medium) loses its 
significance. 

With the potential (2) and the concentrations 
Ca, Cz, the probability of finding a molecule B 
at a distance between ro and ro+dro from a 
molecule A will become proportional to 


Ca Cpe exp ( /kT)ro?dro. (4) 


Adopting as the criterion for the formation of 
X the fact that r should be within these limits 
the rate of reaction becomes** (with k* and K 
constants) : 


k* = KT*(C4Cp/Cz) exp (—¥(r0)/kT). (5) 


Let ko* be the rate constant for some standard 
temperature To, and let the unimportant term 
T'” be neglected in the usual way, then we 
obtain our final formula 


log k* =log ¢1(r0)/kT 
exp (6) 


It does not seem necessary for the present 

purpose to reduce kp to concentration zero, nor 

to a standard value of the dielectric constant. 
From (6) we obtain, for the dependence on 1/T 


R* 
a(t/T) 
exp (7) 
and for the dependence on 1/D 
d log k* 
ap = 
d(1/D) 
= exp (8) 


ar 


= —(1/k)[¢1(r0) 


An inspection of formula (7) will show that a7 
will be negative, whatever the charges 24 and zg 
are, if the non-electrostatic potential (at r=ro) 
is sufficiently large compared with the electro- 
static potential. 

In the case of reactions which conform to 
electrostatics in their dependence on 1/D it 
should be possible, from (7) and (8) and from 
the knowledge of a7 and ap, to derive as well the 
value of the electrostatic potential @z,42s/Dro as 
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TABLE I. Tests on inequalities (9) and (10). 


Reaction Medium i 


ar ap T/D aT/ap 


A. Ions of like sign 


$,0,--+I- MeOH —H,0 313.1 73.12 — 10140 —624 4.27 16.24 

MeOH —H.0 298.1 71.42 — 11000 — 737 4.18 15.24 

B¢B--+OH- EtOH —H.0 298.1 71.42 — 11120 —921 4.18 12.08 
B. Ions of unlike sign 

NH,*++CNO- MeOH —H:20 323.1 63.5 — 9370 260 5.09 — 36.0 

NH,*+CNO- Glycol—H:,0 323.1 63.5 — 10140 184 5.09 —55.2 


TABLE II. A comparison of the non-electrostatic potentials at a distance r=ry with electrostatic pote ntial between ions 


Reaction Medium 


2(r0) X10*8 X10*12 o1(r0) /b2(r0) 


NH,*++CNO- MeOH —H:0 
NH,*+CNO- Glycol —H:O 323.1 
MeOH —H.0 313.1 
B¢B--+OH- EtOH—H:0 298.1 
B¢B--+OH- MeOH —H.0 


63 1.28 1.39 10.9 
73.12 3.66 1.39 3.8 
73.42 5.27 1.52 2.9 


that of the non-electrostatic potential g:(ro) for 
the critical distance. 

Furthermore, it is easy to derive from (7) and 
(8) the following two inequalities which connect 
measurable quantities only: 


_ 1. if A and B have charges of the same sign: 


ar/ap>T/D; (9) 
2. if A and B have charges of opposed sign: 
ar/apn<T/D. (10) 


III. APPLICATIONS 


The first application will be to test the in- 
equalities (9) and (10). Inequality (9) can be 
tested by the data on the persulfate-iodide 
reaction in methyl-water media’ and the brom 
phenol blue-hydroxide reaction in both ethyl 
alcohol-water and methyl alcohol-water media.® 
From Table I it can be seen that a7/ap is three 
to four times as large as T/D. The predictions 
of inequality (10) can be applied to data on the 
ammonium-cyanate reaction in both methyl 
alcohol-water® and glycol-water® media. It is 
evident that in these cases a7/ap is less than 
T/D, to the extent of the inversion of the sign 
of ar/ap. 


7E. S. Amis and J. E. Potts, Jr., J. Am. Chem. Soc. 63, 
2883 (1941). 

8 J. Lander and W. J. Svirbely, J. Am. Chem. Soc. 60, 
1613 (1938). 


In Table II, the potentials g1(ro) and g2(ro) as 
well as their ratios are recorded for different 
sign and valence types of reactants in different 
media. In all cases the rate dependence on 1/T 
was linear as demanded by both the Arrhenius 
equation and by Eq. (7). For bivalent-univalent 
reaction types the value of the non-electrostatic 
potential, ¢:(7o), is roughly three times as great 
as the value of the electrostatic potential, g2(7o), 
while in the uni-univalent type ¢:(ro) is from 
seven to eleven times as great. The electrostatic 
potential would be expected to be larger for the 
higher valence types of reactants. The solvent 
seems to play a rather important part in the 
ratios of these potentials. 

As has been shown above, values of ar and ap 
permit the calculation of gi(7o) and g2(ro) for 
reaction rates obeying electrostatics in their 
dielectric constant dependence. If, however, ap 
is not known, ¢;(7o) may still be calculated from 
ar by use of Eq. (7), provided that ro is known 
from some other source or provided that the 
second term on the right of Eq. (7) is neglected 
as a first approximation. Using this latter method 
of calculation for reactions between ions and 
molecules, molecules and molecules, and reac- 
tants in the gaseous state, the temperature 
dependence of the rates should yield values of 
¢i(ro) substantially like those for reactions be- 
tween ions in solution, since this potential 
arises from repulsive forces between the reactants 
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TABLE III. A comparison of the non-electrostatic potential between reactants calculated from the rate dependence on 
temperature alone and calculated from the rate dependence on both temperature and dielectric constant. 


Reaction 


Value used in 


from rate from rate 
calculating dependence dependence 


column 7 of 


Temp.°C 


B¢B--+OH> 
B¢B--+OH~ 
NH,*++CNO~4 
B¢BOH-+H;0+t 
Sucrose+H;0+t/ 
Sucrose+H;0+/ 


*Fett++4+Sntti 
2N,0/ 
C.H;OCOCL* 
2HI! 

H+H,(P)™ 


EtOH—H,0 
EtOH—H,0 
MeOH 
Glycol—H,0 
MeOH 
EtOH—H.O 
EtOH—H,0 
Dioxane— H,O 
CCl 


Po to ine 


@ See reference 7. 
+ See reference 5. 
¢ See reference 8. 
4 See reference 6. 


¢E. S. Amis and J. B. Price, presented to J. Am. Chem. Soc. for publication. 


/ E. S. Amis and F. C. Holmes, J. Am. Chem. Soc. 63, 2231 (1941). 


9 F, Daniels and E. H. Johnson, J. Am. Chem. Soc. 43, 53 (1921); see F. H. MacDougall, Physical Chemistry (The MacMillan Company, New 


York, 1936), p. 404 


h P. T. Stroup and V. W. Meloche, J. Am. Chem. Soc. 53, 3331 (1931). 


'W. F. Timofeew, G. E. Muchin, and W. G. Gurewitsch, Zeits. f. Physik. Chemie 115, 161 (1925), 

iC. N. Hinshelwood and R. E. Burk, Proc. Roy. Soc. Al06 284 (1924). 

* A. R. Choppin, H. A. Frediani, and G. F. Kirby, Jr., J. Am. Chem. Soc. 61, 3176 (1939). 

'M. Bodenstein, Zeits. f. physik. Chemie 29, 295 (1899); see H. S. Taylor and H. A. Taylor, Elementary Physical Chemistry (D. van Nostrand 


Co., New York, 1937), p. 242. 


™K. H. Geib and P. Harteck, Zeits. f. physik. Chemie, Bodenstein-Festband, 849-62 (1931); see H. S. Taylor and H. A. Taylor, Elementary 


Physical Chemistry (D. van Nostrand Co., New York, 1937), p. 235. 
*k was not extrapolated to ky =0 in these cases. 


at very close distances where the influence of the 
charge and of the dielectric constant is no longer 
appreciable. 

In Table III are recorded the values of ¢;(70) 
calculated in the two ways mentioned above and 
neglecting the second term on the right of Eq. (7) 
in the second method of calculation. The ionic 
examples included in the first method of calcu- 
lation are all the instances recorded in the litera- 
ture where sufficient data is given for this pur- 
pose and where the prerequisite of obedience 
to electrostatic theory is met. It is observable 
from Table III that the values of ¢:(ro) calcu- 
lated from the temperature dependence of the 
rate correspond closely to the values of this 
function calculated from the dependence of the 
rate on both the temperature and dielectric 
constant. Furthermore, the variation between 
the two values is that predicted by electro- 
statics, i.e., g1(70) calculated from ar and ap is 
less than that calculated from ar, neglecting 
the electrostatic term, for ions of like sign, and 


vice versa for ions of unlike sign. It is also evi- 
dent that for all reactions, whether in solution 
or in the gaseous state, and irrespective of 
charge type, the values of ¢1(7o) do not vary by 
an order of magnitude and in many cases agree 
in absolute value. It might be mentioned that, 
in the regeneration by hydrochloric acid of 
brom phenol blue from the colorless carbinol,° 
even though the reaction did not obey the re- 
quirements of electrostatics in its dielectric 
constant dependence of the rate, yet calculation 
of ¢:i(ro) from the rate dependence upon tem- 
perature yielded a value of the non-electrostatic 
potential which compared favorably with that 
for the fading reaction the kinetics of which met 
the requirements of electrostatic theory. 

As to the magnitudes of the non-electrostatic 
potential at the distance of reaction, they are of 
the expected order. Compared with the electro- 
static energies, they are three to ten times as 
large if reasonable values of ro are used. Further- 


9See Table III, reference e. 


|| 
4 Temp. d.c. d.c. x10+2 10412 Media 
8 20-40 55.0-76.7 40 73.12 1.04 
5-45 64.5-78.5 25 71.42 0.93 
5-45 64.5-78.5 25 71.42 1.13 
| 30-60 40.0-63.5 50 63.5 1.70 
; 30-70 35.0-69.9 50 63.5 132 
25-45 
21-41 
21-41 
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more, they can be compared roughly with values 
of the intermolecular energies as calculated from 
the second virial coefficient. Fowler’ has evalu- 
ated the coefficients \ and yp in the well-known 
expression 


gi(r) (11) 


for values of s ranging from 8 to 14 applying 
Eq. (11) to noble gases. Of the reactions given 
above that between the ammonium-ion-cyanate 
ion involves reactants which correspond reason- 
ably well in molecular weight to argon. Using, 


oR, H. Fowler, Statistical Mechanics, second edition 
(Cambridge, 1936), p. 306. 
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then, Fowler’s constants for argon we obtain 
for gi(ro), with 7>=2.5 10-8 cm as determined 
from Eq. (1), values between 1.1 and 3.110-” 
when s ranges from 8 to 14. Thus our value of 
1.7X10-” is of the same order of magnitude. 
The good agreement, however, must be con- 
sidered as fortuitous, as Eq. (11) is very sensi- 
tive to a change in r. For this reason the other 
reactions studied above do not agree nearly as 
well with Fowler’s values, since, for these reac- 
tions, Eq. (1) yields values of ro between 1.2 and 
1.7 X10-*cm which are evidently somewhat small 
if consideration is taken of the large molecular 
weights of the reactants. 
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An X-Ray Investigation of the Structure of Lead Chromate* 
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(Received July 13, 1942) 


An attempt to analyze the structure of lead chromate was undertaken, using the Weissenberg 


x-ray goniometer technique. Positions of the lead and chromium atoms were determined. 
Although the intensities were carefully corrected for absorption, accurate values of the oxygen 
parameters could not be determined from the data at hand. Further work is indicated. There 
is some discussion of the migration of ions in the photoblackening effect which lead chromate 


exhibits. 


RILL! and others*:* have established that 

lead chromate (crocoite) is monoclinic with 
cell dimensions a=7.10A, b=7.40A, c=6.80A, 
8=102°27’. There are 4 molecules (PbCrO,) per 


cell. The probable space group is Con—P2,/n. 
These data do not permit lead chromate to be 
classified with any of the previously determined 
structures in the oxyacid group of substances,‘ so 
that a detailed x-ray investigation would seem 
to be of interest. In addition, a knowledge of the 
structure is important for an interpretation of 


* Submitted to the faculty of Bryn Mawr College in 
partial fulfillment of the requirements for the degree of 
doctor of philosophy. Gratitude is expressed to Dr. A. L. 
Patterson for his guidance and direction of the investiga- 
tion. A grant made to Bryn Mawr College by the Elizabeth 
Thompson Science Fund made possible the purchase of 
much of the apparatus used. 

1R. Brill, Zeits. f. Krist. 77, 506 (1931). 

? Gossner and Mussgnug, Zeits. f. Krist. 75, 410 (1939). 

3S. von Glisczcynski, Zeits. f. Krist. 101, 1 (1939). 

4See R. W. res Wyckoff, The Structure of Crystals 
(Chemical Catalogue Company, New York, 1931), second 
edition, Chap. 16. 


the photoblackening process which lead chromate 
exhibits. 

Simple rotation and Weissenberg equi-inclina- 
tion’ x-ray photographs were taken about the c 
and b axes. The indexing of these photographs 
verified the lattice and space group previously 
established. Unfiltered copper radiation was used 
throughout. Intensities were visually estimated. 
To aid in putting the intensities on a numerical 


TABLE I. 
u v w 
(fractions of cell 
Distance vector Projection translations) 

(2x 2y 22) c face 0.43 0.30 

b face 0.45 0.80 
(3 3—2y 3) c face 3 0.20 

b face 3 
(3—2x 3 3—2z) c face 0.07 

b face 0.05 0.70 


5M. J. Buerger, Zeits. f. Krist. 88, 356 (1934). 
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STRUCTURE OF LEAD CHROMATE 


scale, two exposures of different times were made TABLE III. Observed and calculated structure factors. 


in 
ed at each setting. Use was also made of the relative 
“12 intensity of reflections from the Cu Ka and Kp 
of components. 200 360 —340 200 390 — 343 
le. The extremely high absorption coefficient of = 
n- lead chromate (u=1215/cm for the radiation 800 0 - § 800 0 2 
si- used) makes it necessary to be very careful in 
er correcting the intensities for absorption. In pre- 
1 1 iil 1 
as liminary work, crystals were used in the form 519 175 133 501 160 163 
ic- in which they occur in nature, and plane by 610 200 — 183 501 0 — 45 
nd plane absorption corrections based on the = 
all Hendershot® method were applied. Even with 020 73 —117 901 30 — 106 
lar refinements in the method of making the cor- = = 
rection, inconsistencies in the results from these 320 280 246 302 290 —224 
crystals indicated that the best natural Is 420 0 — 61 402 130 —127 
which cou ound’ lac the perfection of 60 0 27 602 105 173 
form required for a safe application of this = 
method. The artificial crystals available were too 139 46 -~ 63 302 130 ~ 201 
small to give satisfactory photographs. It was = - 135 
therefore decided to prepare cylindrical speci- 439 “59 61 303 110 116 
mens. 530 210 —201 303 175 — 182 
A few good single crystals of crocoite of suf- 9 34 a 50380 == 
730 165 202 503 +120 135 
ficient size were selected. The crystal to be 830 22 62 = 20 3 
H H 040 8 243 — 248 03 80 — 83 
ground was mounted ina goniometer head and 33 
set along the appropriate axis with the aid of a 249 240 017 004 «250 —242 
two-circle goniometer. It was then transferred ree 
with little disturbance of the orientation to a 549 35 as 404 135 ~232 
rigid mount in the headstock of a small lathe. 640 30 66 404 90 - 47 
A small piece of sandpaper attached at one end ~ 
ate to a tool blank and free at the other end served 150 0 0 804 60 116 
as a grinding tool. A high point on the sandpaper 
1a- was the actual cutting tool. This paper is stiff 450 155 —172 305-35 — 69 
ec enough to furnish the pressure required to cut 1 - = 
yhs the soft but brittle crystal, and flexible enough to 759 0 0 505 ~=—s- 80 117 
sly prevent the frequent breakage attendant upon 
sed the use of a rigid tool. The finished cylinders 69 295 183 206 40 112 
ed. were about mm in diameter and 1-2 mm long. = = 
cal Intensities of reflections from the cylindrical 569 53 68 406 215 230 
TABLE II, Parameters in fractions of cell translations. pos 
107 
w Ps y s 370 65 —118 307 60 —129 
307 
= Pb 0.220 0.150 
Cr .215 .150 
).80 
All values +0.004. 
2 °Q. P. Hendershot, Rev. Sci. Inst. 8, 324 (1937). A dis- 80 31 43 
cussion of the application of this method will appear 4 
shortly. 190 0 — 21 
).70 290 25 — 63 


’ Crocolte from Tasmania, Brazil, and the Urals was 

ordon, an r. E. atson. Artificial crystals were Paramet sed in calculating: * Pb =(0.219, 0.150), C 
obtained from Dr. S. é Horning. =(0.216.0.150); **Pb =(. 220, (x2, £2) 
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crystals, corrected for absorption by the method 
of Bradley® were used in the main part of the 
work. 

The general position of P2,/n is (x y z)(# 9 2) 
metry of the F? series is P2/m. A set of atoms of 
atomic number Z in such a general position will 
produce an interatomic distance peak in the F? 
series® of magnitude Z? at each of the general 
points 4: (0) 2x 2y 22 of the space group P2/m,!° 
and a peak of magnitude 2Z? at each of the special 
points 2: (1) } 3+2y 3 and 2: (nm) $+2x 4 4+22 of 
the same group. If two general points x1y1z; and 
X2yo%2 Of P2,/n are occupied, peaks of equal mag- 
nitude are produced at the general points 4: 
(0) S,S,Sz, 4: (0) D:D,D., 4: (0) +S, 
4: (0) $+D, 4+S, 4+D. of P2/m, where 
S,=Xet%1, D,=%X2—-%1, etc. 

Using the intensities of the (4kO) and the 
(hOl) reflections, properly corrected for polar- 
ization and absorption, the projections of the F? 
series on the b face and c face were computed and 
plotted. From these maps and from symmetry 
considerations it was concluded that all atoms 
must occupy general positions. The only reason- 
able interpretation of the strongest peaks, cor- 
responding to the lead-lead distances, is given in 
Table I. This gives a unique set of parameters 
for the lead atoms, x,;=0.22, y;=0.15, 2:=0.40. 

Because of the high atomic number of lead it 
can be assumed that the signs of the structure 
factors are in the main determined by the con- 
tributions of the lead atoms, Fp,(hk/). The signs 
of the contributions Fp,(hk0O) and Fp,(h0l) were 
accordingly determined from the above lead 
parameters. These signs were combined with the 
experimentally determined numerical values of 
the total structure factors and the set of trial 
F’s thus constructed was used in synthesizing 
the Fourier series for the projected density on 
the ¢ and 6 faces. The chromium atoms were 
readily located on these maps. Parameters of the 
lead and chromium atoms obtained from the 
electron density maps are listed in Table II. The 
position of the lead atoms is identical with that 


8 A. J. Bradley, Proc. Phys. Soc. 47, 879 (1935). 


® A. L. Patterson, Zeits. f. Krist. 90, 517 (1935). 
1 Notation of Int. Taf. sur Best. der Krist. 
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obtained from the F? maps. All peaks due to the 
lead-chromium distances appear in the right 
positions and orders of magnitude on the F? 
series maps. The structure factors calculated 
from these lead and chromium parameters are in 
good agreement with the observed magnitudes 
of the F’s (Table III). No sign changes are 
introduced by the contributions of the chromium 
atoms. Thus the assumption on which the calcu- 
lation of the F series was based is justified. 

No values are listed for the oxygen parameters 
since they cannot be located accurately from the 
data at hand. The maps give some indication of 
the positions of the oxygen atoms but the 
presence of spurious peaks on the maps, probably 
due to cutting off the Fourier series too soon, 
indicates that caution must be used in inter- 
preting the weaker peaks. An extension of the 
data to higher order reflections would be desir- 
able. Although plausible oxygen parameters can 
be obtained by adjusting the calculated structure 
factors, this procedure seems worthless without 
the support of direct evidence since the structure 
factors are very insensitive to large variations in 
the oxygen parameters. With the aid of a scale 
model attempts were made to find possible struc- 
tures by packing nearly tetrahedral chromate 
groups about the lead ions. It is difficult to find 
adequate radii for the lead and chromate radicals 
in the literature and these attempts were incon- 
clusive. From the lead and chromium positions 
and from indications of possible oxygen arrange- 
ments it seems that the structure is fairly close 
packed. The coordination number of oxygen 
atoms about lead is probably 8 or more. 

In the light of these findings it would appear 
that the mechanism of photo-conductivity and 
photoblackening in lead chromate cannot depend 
on the migration of ions through a perfect lattice. 
Migration by means of holes due to lattice defects 
remains as a possibility. If any migration in a 
perfect lattice takes place it probably is accom- 
panied by a breakdown of chromate groups and 
therefore of the lattice as a whole." 


11T am indebted to Dr. Frederick Seitz for a discussion 
of this interpretation. It was Dr. Seitz’ interest in the 
problem which suggested the research reported here. 
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Letters to the Editor 


HIS section will accept reports of new work, provided 
these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 


The Carbon-Carbon Bond Strengths in Ethane, 
Propane, and n-Butane 


E. R. VAN ARTSDALEN 
Department of Chemistry, Lafayette College, Easton, Pennsylvania 
August 19, 1942 


ECENTLY the carbon-hydrogen bond strengths in 
methane and ethane were determined from kinetic 
data by Andersen, Kistiakowsky, and Van Artsdalen' who 
gave values which agree very well with those obtained by 
Stevenson? from measurements of electron impact po- 
tentials. In the calculation of the carbon-hydrogen bond 
strengths the various energy relations were not corrected to 
absolute zero. When this is done, one finds that the CH;— 
bond strength, in round numbers, is 101 kcal. at 0°K, 
102 kcal. at room temperature, as reported, and 103 kcal. 
at 450-475°K. 

The heat of dissociation of carbon-hydrogen bonds may 
be combined with other established thermal data to cal- 
culate carbon-carbon bond strengths in several simple 
hydrocarbons. From thermal data tabulated by Rossini® 
one obtains the following at 298.16°K: 


2CH, = C:H.+H2—- 15.54 kcal. 


This can be combined with the value of 102+1 kcal. for 
the strength of the CH;—H bond and with the value for 
the heat of dissociation of hydrogen‘ to obtain a value 
for the CH;— CH; bond strength. Thus at 298.16°K 


2CH,=C2He +H, — 15.5 kcal. 
2CH,=2CH;+2H—204 kcal. 
H.=2H — 102.9 kcal. 


C.H,g=2CH; — 85.6 kcal. 


The probable error is +2.3 kcal. 

This value is in good agreement with that calculated 
by Stevenson.? Trenner, Morikawa, and Taylor’ obtained 
97.6 kcal. employing their value of 108 kcal. for the 
CH;—H bond strength. It seems certain that the carbon- 
carbon bond strength in ethane must be lower than the 
carbon-hydrogen bond strength in ethane (98+2 kcal.) to 
account for pyrolytic products. Rice and Johnston’ esti- 
mated from activation energies of free radical formation 
at high temperature a C—C bond of 74 kcal. Because 
chain decomposition appears to be less frequent than had 
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been thought, Kistiakowsky’ suggested that the true bond 
strength is higher than 74 kcal. 

Similar calculations may be made to ree the 
carbon-carbon bond strength in propane and the center 
bond in normal butane. The necessary data (heats of 
formation of the hydrocarbons, C—H bond strengths in 
CH, and CsH¢, and the heat of dissociation of He) are 
given in the above references. This method of calculating 
carbon-carbon bonds may be applied wherever the appro- 
priate carbon-hydrogen bonds and heats of formation of 
the hydrocarbons are known. The following resumé of 
values for 298.16°K is obtained: 


CH; —CH; 85.6+2.3 kcal. 
C.H;—CH; 83.8+3.3 kcal. 
C2H;—CoH; 82.4+4.5 kcal. 


Lowering of the C—C bond strength with increasing 
chain length is observed in agreement with expectation. 
It is noted that a CH: increment lowers the bond strength 
about 1.5 kcal. However, qualitative rather than quanti- 
tative significance should be attached to this because of 
the rather large probable errors. 


usin Kistiakowsky, and Van Artsdalen, J. Chem. Phys. 10, 
194 

2D. P. Stevenson, J. Chem. Phys. p. 291 (1942). 

3F. D. Rossini, Chem. Rev. 27, 1 40). 

4 Bichowsky and Rossini, (Reinhold, 1936). 

6 Trenner, Morikawa, and Phys. 5, 203 (1937). 

s Rice and Johnston, z. Am. Chem. Soc. 56, 214 (1934). 

7G. B. Kistiakowsky, J. Phys. Chem. 41, 175 (1937). 


The Partition Function of a Gas of 
Hard Elastic Spheres 


O. K. RIcE 
University of North Carolina, Chapel Hill, North Carolina 
September 8, 1942 


ONKS'! has derived the equation of state of a one- 
dimensional gas consisting of N hard elastic spheres 

of radius ¢, constrained to move along a line of length /. 
The partition function, found readily from his work, is 


(1) 


where No/I. The factor occurring because the whole 
effective length 1(1—6) can be shared, gives rise to a 
“communal entropy” of Nk. 

In three dimensions Tonks considers separately the 
dilute and the concentrated gas. For the former he uses 
an equation of state of the virial type, Eq. T(31). (T 
refers to Tonks’ equations.) For the concentrated gas, 
Tonks obtains an equation of state, Eq. T(32), by an’ 
ingenious argument, which is equivalent to using as the 
partition function 


— 04) PY, - (2) 


essentially the cube of Eq. (1), with v! (v=volume per 
molecule) taking the place of 1/N, and 6! replacing @. In 
Eq. (2), thus the free volume v(1— 6!) vanishes 
when the spheres are in contact in a face-centered cubic 
lattice. To bridge the gap between Eq. T(31) and Eq. 
T(32), Tonks gives an interpolation formula, Eq. T(33). 
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In Eq. (2) the factor e*” corresponds to the large com- 
munal entropy 3kN, because lengths are shared in the 
three directions of space, rather than volume being shared 
as such. That such a large communal entropy can exist 
does not seem to have been explicitly pointed out before, 
though it is implicit in an earlier note by the writer.? In 
the dilute gas, the communal entropy is only kN. Thus 
the change of entropy per mole from a small volume 1 
to a large volume 1 will be 


AS=RT In (3) 


AS can also be calculated by integrating the thermo- 
dynamic equation v(0p/0T),=(0S/dInv)7, using Eq. 
T(33) (see Fig. 1). This gives a AS about 1 entropy unit 
smaller than does Eq. (3), despite the —2R in the latter. 
A better interpolation is shown in Fig. 1; this assumes that 
Eq. (2) holds to somewhat lower values of 6 than does 
use of Eq. T(33). The agreement that can thus be obtained 
shows the essential consistency of Tonks’ ideas and seems 
to offer support to the correctness of the large communal 
entropy inherent in Eq. (2). 


iL. Tonks, Phys. Rev. 50, 955 (1936 
20. K. Rice, J. Chem. Phys. 9, 121 M041). 


The Partition Function of a Simple Liquid 


O. K. RICE 
University of North Carolina, Chapel Hill, North Carolina 
September 8, 1942 


F we neglect the mutual potential energy of the atoms 
and consider only the thermal energy, a face-centered 
cubic solid bears a resemblance to a close-packed gas of 
hard spheres. We can get an effective value 69 for @ of 
Eq. (2) (preceding note) by equating Q to the thermal 
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part of the partition function for the solid, e7?/@* (where 
© is the Debye characteristic temperature). @ may be 
obtained as a function of interatomic distance, if the inter- 
molecular potential energy is known. 

This may be extended to the liquid, if we consider it to 
be formed from the solid by introduction of holes through 
removal of atoms; these holes then contribute to the free 
volume.! If there are N atoms and 2 holes, we write 
0=0,.N/(n+N), where 9 is the value for the solid with 
the same interatomic distance, obtained as in the pre- 
ceding paragraph. We substitute @ into Eq. (2), noting 
also that the molecular volume has been increased by a 
factor (n+N)/N. 

This procedure gives only that part of the partition 
function due to the thermal motion. We must consider, also, 
the energy necessary in order to produce a hole, already 
discussed in previous work.? But, in addition, in order 
that these holes should actually be used, it is necessary 
for the atoms to move past each other, with fluctuations 
in the distances between nearest neighbors. These fluctua- 
tions show up in the atomic distribution curves obtained 
from x-ray investigations, which show a broader band of 
nearest-neighbor distances than can be accounted for by 
thermal motion alone. This results in a further increase 
in the average potential energy, whose order of magnitude 
in the liquid can be estimated from the atomic distribution 
curves. In order that the holes should be effectively used, 
the increase in entropy gained therefrom must more than 
counterbalance the increase of energy due to both causes 
noted above. 

These ideas have been applied to argon. In doing this 
the free energy is minimized using n/N as an additional 
parameter. The picture appears to work very well in ex- 
plaining the properties of the liquid at the normal melting 
point. At 193°K, the theory makes the liquid appear to be 
too compressible. This seems to be due in large part to 
two causes. (1) The potential energy curve previously 
deduced® and used in these calculations probably does 
not rise steeply enough at small distances. (2) It will be 
recognized that if a close-packed model gas of hard elastic 
spheres contains only a few holes, these holes cannot be 
used to increase the freedom of motion, since they are 
blocked off by surrounding spheres. As the number of 
holes increases, a point will rather suddenly be reached 
where they can be shared. In our model of a liquid, some 
similar effect will occur. This matter will be discussed in 
a fuller report, in the light of the data and the detailed 
calculations. 

The entropy of the liquid obtained on this model agrees 
fairly well with that calculated by the method given in a 
previous note.! It is slightly smaller, which is perhaps in 
the wrong direction, but it increases somewhat more slowly 
with 2/N, and in this respect appears to give a better 
indication of the behavior of the liquid. 

1 Compare O. K. Rice, J. Chem. Phys. 5. St (1941). 


20. K. Rice, J. Chem. Phys. 7, 883 (19 
30. K. Rice, J. Am. Chem. Soc. 63, 3 (194i). 
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